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Abstract 

We study the heat capacity and neutrino emission reactions (direct and modified Urea processes, 
nucleon-nucleon bremsstrahlung, Cooper pairing of nucleons) in matter of supranuclear density 
of the neutron star cores with superfluid neutrons and protons. Various superfluidity types are 
analysed (singlet-state pairing and two types of triplet-state pairing, without and with nodes of 
the gap at a nucleon Fermi surface). The results are used for cooling simulations of isolated 
neutron stars. Both, the standard cooling and the cooling enhanced by the direct Urea process, 
are strongly affected by nucleon superfluidity. Comparison of cooling theory of isolated neutron 
stars with observations of their thermal radiation may give stringent constraints on the critical 
temperatures of the neutron and proton superfluidities in the neutron star cores. 
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1 Introduction 

Neutron stars (NSs) are unique astrophysical objects. First, their observational manifestations are 
numerous (radio and X-ray pulsars, X-ray bursters, X-ray transients, etc.; Shapiro and Teukolsky, 
1983). Second, their structure and evolution are determined by properties of matter under extreme 
conditions which cannot be reproduced in laboratory (supranuclear densities, superstrong magnetic 
fields, superfluidity of baryon component of superdense matter, etc.). 

It is widely accepted (e.g., Imshennik and Nadyozhin, 1982; Shapiro and Teukolsky, 1983) that NSs 
are born at the final stage of evolution of normal stars of mass M > 8Mq in gravitational collapse of 
their cores {Mq is the solar mass). During collapse, matter of central layers is compressed to nuclear 
densities and is enriched by neutrons. As a result, a compact NS is created of mass M ~ Mq and 
radius 7?. ~ 10 km. Its core consists mainly of neutrons (with some admixture of protons, electrons, 
and - possibly - hyperons and other particles). A NS is born hot, with the internal temperature about 
10^^ K, but cools down rapidly due to powerful neutrino emission from the internal layers. 

NSs possess strong gravitational fields; gravitational acceleration on their surfaces is as large as 
~ (2-3) X 10^^ cm s~^. The stellar radius TZ is only 2-3 times larger than the gravitational radius 
TZg = 2GM /(? ^ where G is the gravitational constant and c is the light speed. Therefore, the effects 
of General Relativity are important in NS life. 

In this review, we restrict ourselves by consideration of structure and thermal evolution of isolated 
NSs of age t < 10^ yr. Evolution of NSs in close binaries is more complicated; it is described, for 
instance, by Lewin et al. (1995). 

Superfluidity and superconductivity of nucleons in NS interiors play important role in evolution of 
isolated NSs. For instance, sudden changes of spin periods (glitches) demonstrated by some pulsars 
are commonly explained by interaction of normal and superfluid components of matter in the stellar 
crusts (Pines, 1991). Superfluidity affects the heat capacity and neutrino luminosity of the star, and 
therefore its cooling. The effect of superfluidity on the cooling is the main subject of the present 
review. 

At a certain cooling stage (10^-10^ yr) the NS internal temperature can strongly depend on critical 
temperatures of transitions of nucleons into superfluid state. This feature was first mentioned by Page 
and Applegate (1992) who called NSs "thermometers" for measuring the nucleon critical temperatures 
in asymmetric nuclear matter. Microscopic calculations of such temperatures are complicated owing to 
the absence of explicit manybody relativistic quantum theory which would describe adequately strong 
interactions of particles of various species. However, the critical temperatures can be studied by the 
astrophysical method, by confronting cooling simulations with observations of thermal radiation of 
isolated NSs. This work is not finished yet but certain results are already obtained. The aim of this 
review is to describe the above method and main results obtained to-date (by January 1999). 



2 Overall properties of neutron stars 

2.1 Structure 

If standard measurements are used, neutron stars (NSs) are hot objects. In several hundred thousand 
years after their birth their internal temperature exceeds 10^ K. However, the equation of state of 
internal NS layers is practically independent of temperature since the main contribution into pressure 
comes from strongly degenerate fermions of high energy (Shapiro and Teukolsky, 1983). 

A NS can be subdivided into the atmosphere and four internal regions: the outer crust, the inner 
crust, the outer core, and the inner core as shown in Fig. |^. 

The atmosphere is a thin plasma layer, where the spectrum of thermal electromagnetic NS 
radiation is formed. In principle, this radiation contains valuable information on stellar parameters 
(on temperature, gravitational acceleration and chemical composition of the surface, on magnetic field, 
etc., see Sect. 8.1) and, as a result, on the internal structure. Geometrical depth of the atmosphere 
varies from some ten centimeters in a hot NS down to some millimeters in a cold one. Very cold NSs 
may have no atmosphere at all but a solid surface. 

NS atmospheres have been studied theoretically by many authors (see, e.g., review articles by 
Pavlov et al., 1995; Pavlov and Zavlin, 1998, and references therein). Construction of the atmosphere 
models, especially for the cold NSs (with the surface temperature T^ < 10^ K) with strong magnetic 
fields 10^^-10^^ G, is far from being complete owing to complexity of calculations of the equation of 
state and spectral opacity of the atmospheric plasma (Sect. 8.1.2). 

The outer core (outer envelope) extends from the atmosphere bottom to the layer of the 
density p = Pd ^ 4.3 x 10^^ g cm^^ and has a depth about some hundred meters (Shapiro and 
Teukolsky 1983). Its matter consists of ions and electrons. A very thin (no more than several meters 
in a hot NS) surface layer contains a non-degenerate electron gas. In deeper layers, electrons constitute 
a strongly degenerate, almost ideal gas, which becomes relativistic at p ':$> 10^ g cm~^. For p > 10^ 
g cm~^, atoms are fully ionized by the electron pressure, being actually bare atomic nuclei. The 
electron Fermi energy grows up with increasing p, the nuclei suffer beta-decays and are enriched by 
neutrons (see, e.g., Haensel and Pichon, 1994). At the base of the outer core (p = pd) neutrons start 
to drip from nuclei producing a free-neutron gas. 

The depth of the inner crust (inner envelope) may be as large as several kilometers. The 
density p in the inner crust varies from pd at the upper boundary to ^ 0.5po ^^ the base. Here, 
po = 2.8 X 10^^ g cm~^ is the saturation nuclear matter density. Matter of the inner crust consists 
of electrons, free neutrons and neutron-rich atomic nuclei (Negele and Vautherin, 1973; Pethick and 
Ravenhall, 1995). The fraction of free neutrons increases with growing p. At the core bottom (in 
the density range from 10^^ to about 1.5 x 10^^ g cm~^) the nuclei may form clusters and have non- 
spherical shapes (Lorenz et al., 1993; Pethick and Ravenhall, 1995). The nuclei disappear at the 
crust-core interface. Neutrons in the inner crust may be in a superfluid state (Sect. 3.1). 

The outer core occupies the density range 0.5/9o ^ P ^ "^Po ^md can be several kilometers in 
depth. Its matter consists of neutrons n with some (several per cent by particle number) admixture 
of protons p and electrons e (the so called standard nuclear composition). Composition of matter is 
determined by the conditions of electric neutrality and beta-equilibrium with respect to the reactions 
n— >p + e-|-Pe; p + e ^>- n + Ve, where Ve and i>e stand for electron neutrino and antineutrino, 
respectively. The electric neutrality requires equality of the electron and proton number densities: 
Tip = rig. The beta equilibrium establishes a relationship between chemical potentials of the particles: 
Pn = Pp + Pe- The neutrino chemical potential is omitted here since immediately after the birth a 
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Figure 1: Schematic cross section of a 1.4 M© neutron star. Stellar parameters depend strongly on the equation of 
state of internal layers. 



NS becomes fully transparent for neutrinos which freely escape the star (Friman and Maxwell, 1979; 
Shapiro and Teukolsky, 1983). 

All npe-plasma components are strongly degenerate. The electrons form almost ideal relativistic 
Fermi gas. The neutrons and protons, which interact via nuclear forces, constitute a strongly non-ideal 
non-relativistic Fermi-liquid. It can be described using the formalism of quasi-particles (Botermans 
and Malfliet, 1990; Baym and Pethick, 1991). The Fermi momentum ppj of particle species j is related 
to their number density Uj as ppj = h (Svr^ Uj)^'^. Therefore, the electric neutrality implies the equality 
of the Fermi momenta of protons and electrons: pp = p-p^. 

Calculations show (e.g., Shapiro and Teukolsky, 1983) that for p ~ po the chemical potentials 
of neutrons and electrons are p„ ~ /Xe ~ 60-100 MeV, while the chemical potential of protons is 
much lower, jXp ~ 3~6 MeV. With increasing density, the Fermi energies of particles grow up, so that 
new particles can be created. First of all, these are muons which, like electrons, constitute a strongly 
degenerate, almost ideal gas. The appearance of particles in the inner NS core will be discussed below. 

Let us emphasize that properties of the NS crust (p < 0.5po) a-re described by sufficiently reliable 
microscopic theories. The situation with matter of supranuclear density, p > po, is quite different. 
Laboratory data on properties of this matter are incomplete, and reliability of the theories decreases 
with growing p. Exact selfconsistent quantum theory of matter of supranuclear density has not been 
constructed yet. Many theoretical equations of state have been proposed which can be subdivided 
conventionally into the soft, moderate and stiff ones with respect to compressibility of matter. These 
equations of state are considerably different for p > po . For instance, at p « 4po a soft equation 
of state proposed by Pandharipande (1971) (based on the Reid, 1968, nucleon-nucleon soft-core 
potential) gives the pressure which is about one order of magnitude lower than that given by a stiff 
equation of state constructed by Pandharipande and Smith (1975) using the mean field approach. 

Almost all microscopic theories predict the appearance of neutron and proton superfluidity in the 



outer NS core (Sect. 3.1). The proton superfluidity is accompanied by superconductivity (most likely, 
of second type) and affects evolution of the internal magnetic fields (see, e.g., Ruderman, 1991). 

The outer core of a low-massive NS extends to the stellar center. More massive NSs possess also 
the inner core. Its radius may reach several kilometers, and its central density may be as high 
as (10-15)po- Composition and equation of state of the inner core are poorly known. It should be 
emphasize that these properties are most crucial for determining structure and evolution of massive 
NSs, and they constitute the main NS "mystery" . Several hypotheses have been put forward and are 
being discussed in the literature; it is impossible to reject any of them at present: 

(1) Hyperonization of matter — appearance of S- and A-hyperons (see, e.g., Shapiro and Teukolsky, 
1983). The fractions of p and e may be so high that the powerful direct Urea process of neutrino emis- 
sion becomes allowed (Lattimer et al., 1991) as well as similar reactions involving hyperons (Prakash 
et al., 1992). In these cases, the neutrino luminosity of the star is enhanced by 5-6 orders of magnitude 
(see, e.g., Pethick, 1992) as compared to the standard neutrino luminosity of the outer core produced 
mainly by the modified Urea process. This accelerates considerably NS cooling (Sects. 2.2 and 7). 

(2) The second hypothesis, proposed in different forms by Bahcall and Wolf (1965a, b), Migdal 
(1971), Sawyer (1972), and Scalapino (1972), assumes the appearance of pion condensation. The 
condensates soften the equation of state and enhance the neutrino luminosity by allowing the reactions 
of the direct Urea type (Bahcall and Wolf, 1965a, b; Maxwell et al., 1977; Muto and Tatsumi, 1988). 
However, many modern microscopic theories of dense matter predict weakly polarized pionic degrees 
of freedom which are not in favor of pion condensation (Pethick, 1992). 

(3) The third hypothesis predicts a phase transition to the strange quark matter composed of 
almost free u, d and s quarks with small admixture of electrons (see, e.g., Haensel, 1987; Schaab et al., 
1996). In these cases the neutrino luminosity is thought to be considerably higher than the standard 
luminosity due to switching on the direct Urea processes involving quarks (Iwamoto, 1980; Burrows, 
1980; Pethick, 1992). However, in some models (e.g., Schaab et al., 1997a) the presence of quarks does 
not enhance the neutrino luminosity. 

(4) Following Kaplan and Nelson (1986), Nelson and Kaplan (1987) and Brown et al. (1988), 
several authors considered the hypothesis on kaon condensation in dense matter. Kaon condensates 
may also enhance the neutrino luminosity by several orders of magnitude (Pethick, 1992). A critical 
analysis of contemporary theories of kaon condensation was given by Pandharipande et al. (1995). 

For any equation of state, one can build a set of NS models with different central densities pc- 
This can be done by solving numerically the equation of hydrostatic equilibrium with account for the 
effects of General Relativity (the Oppenheimer-Volkoff equation, see, e.g., Shapiro and Teukolsky, 
1983). With increasing p^-, the stellar mass M usually grows up, and the radius TZ decreases (a star 
becomes more compact). As a rule, the mass M reaches maximum at some po which corresponds to 
the most compact stable stellar configuration. The mass M^ of this configuration is the maximum NS 
mass, for a given equation of state. Stellar models with higher central density are usually unstable 
(with respect to collapsing into a black hole) and cannot exist for a long time. Using a set of NS 
models, one can construct a "mass-radius" diagram which depends strongly on the equation of state 
in the NS core. The softer the equation of state, the more compact NS, and the lower the maximum 
mass. For instance, this mass ranges from ~1.4 to ~1.6 Mq for different soft equations of state; from 
~1.6 to ~1.8 Mq for moderate equations of state, and from ~1.8 to ~3 Mq for stiff equations of 
state. The problem of maximum NS mass is crucial for identification of black holes in binary systems 
(Cherepashchuk, 1996). 

The equation of state in the NS cores can be studied by confronting theory with observations 
in different ways. The majority of methods are based on determination (constraint) of NS mass 



and/or radius and comparison of these observational results with the mass-radius diagrams for different 
equations of state (see, e.g., Shapiro and Teukolsky, 1983). Unfortunately, no absolutely decisive 
arguments have been given so far in favor of the stiff, moderate or soft equations of state. One can 
definitely rule out only the ultra-soft equations of state which give the maximum NS mass lower than 
1.44 Mq, the mass of the Hulse-Taylor pulsar (PSR 1913+16), which is the most massive known star in 
close double neutron-star binaries (where NS masses are determined very accurately). In this review, 
we will discuss another method to explore properties of superdense matter — by NS cooling. 

2.2 Neutrino emission 

In ~ 20 s after its birth, a NS becomes transparent for neutrinos (Shapiro and Teukolsky, 1983) 
generated in its interiors. Leaving the star, neutrinos carry away energy and cool the star. Therefore, 
a study of neutrino reactions is important for the NS cooling theories. The most powerful neutrino 
emission is produced in the NS core. Thus we will not discuss the neutrino reactions in the NS 
crust. Their detailed description can be found, for instance, in the review article by Imshennik and 
Nadyozhin (1988) as well as in the articles by Soyer and Brown (1979), Itoh et al. (1984, 1996), Pethick 
and Thorsson (1994), and Kaminker et al. (1997, 1999a). Typical neutrino energies are much higher 
than their assumed rest mass energies; therefore, neutrinos can be considered usually as massless 
particles. 

There are many neutrino reactions in the NS core (see, e.g., Shapiro and Teukolsky, 1983; Pethick, 
1992). Their efficiency is determined by the equation of state of matter which is poorly known (see 
above). For certainty, we will not consider exotic models of strange matter, pion or kaon condensates, 
but outline the main neutrino reactions in npe matter. A more detailed analysis of these reactions, 
particularly in the presence of nucleon superfiuidity, is given in Sects. 4-6. 

Modified Urea proeess. Starting from the classical article by Chiu and Salpeter (1964) this 
processes has been treated as the main neutrino generation process for the standard NS cooling. 
Cooling is called standard provided very powerful neutrino reactions are absent, like the direct Urea 
process or analogous processes in exotic models of dense matter. 

The modified Urea reaction is similar to the familiar reactions of beta decay and beta capture, but 
involves additional spectacular nucleon. In the npe-matter, the reaction can go through two channels 

n + n^n + p + e + i'e, n+p + e^n + n + Ue; (1) 

n + p^p + p + e + Ve, p + p + e^n + p + Ue, (2) 

which we define as the neutron and proton branches of the modified Urea process, respectively. 

The spectacular nucleon is needed to conserve momentum of reacting particles. The familiar 
beta-decay and beta-capture reactions (called the direct Urea process discussed below) are forbidden 
(strongly suppressed) in the outer NS core. The suppression is due to relatively low fractions of elec- 
trons and protons, owing to which the Fermi momenta of n, p, and e do not obey the "triangle rule", 
PPn — PFe~^PFp' required for momentum conservation (the momenta of emitting neutrinos Pu ^ k^T/c 
can be neglected: they are determined by temperature of matter, T, and are much smaller than the 
momenta of other particles; here k-Q is the Boltzmann constant). The neutrino energy emission rate 
(emissivity) Q^ "Mn the neutron branch ([l|) of the modified Urea process for p ~ po can be estimated 
as (Friman and Maxwell, 1979) Q^^""^ ~ 10^2 Tl erg cm'^ s-^ where Tg = T/IO^ K. The proton 
branch is nearly as efficient: Q^^^^ ~ Q^^") (Sect. 5.1). 



Neutrino bremsstrahlung due to nucleon-nucleon scattering. The standard neutrino lu- 
minosity is determined also by the processes of neutrino bremsstrahlung radiation in nucleon-nucleon 
collisions 

n + n^n + n + v + r', n+p^n+p + v + v, p + p^p + p + v + v, (3) 

allowed in the entire NS core. Here v means neutrino of any flavor (fg, f^, or Ut). In normal (non- 
superfluid) matter the bremsstrahlung processes (Friman and Maxwell, 1979) are weaker than the 
modified Urea process: Q^^^' ~ (10^^-10^'^) Tg erg cm~^ s~^. However, they can be important in 
superfluid matter (Sect. 5). 

Direct Urea process. As mentioned above, a sequence of beta-decay and beta-capture reactions, 

n^p + e + i>e-,P + e.^n + Ve., (4) 

called the direct Urea process, is forbidden in the outer NS core due to insufficiently high fractions of 
e and p. It has been thought for a long time that the process is forbidden also in the inner NS core. 

However, the process (Q) becomes allowed (Lattimer et al., 1991), if the fraction of protons (among 
all baryons) Xp = Up/rnj exceeds some critical value Xp = Xc- In npe-matter, this happens for pp^ < 
2pj^, which gives Xc = 1/9 = 0.111. If muons are present for the same number density of baryons 
nf), the proton fraction appears to be slightly higher than in the npe matter, and the electron fraction 
slightly lower. In this case, x^ is higher and reaches 0.148 for relativistic muons (Lattimer et al., 1991). 

In the simplest model of superdense matter as a gas of non-interacting Fermi particles (Shapiro 
and Teukolsky, 1983) the proton fraction is not high: Xp < Xc for any density. However, this may be 
not so for realistic equations of state. This circumstance was first outlined by Boguta (1981) but his 
article has been unnoticed for a long time. 

It was the paper by Lattimer et al. (1991) which initiated wide discussion of the direct Urea process. 
The authors showed that, for many realistic models of matter, Xp exceeded slightly Xc at densities 
several times higher than the standard nuclear matter density. Therefore, the nucleon direct Urea 
process can be allowed in the inner cores of rather massive NSs. Moreover, practically all equations 
of state in the inner stellar core, which predict appearance of hyperons, open direct Urea reactions 
involving hyperons (Prakash et al., 1992). 

According to Lattimer et al. (1991) the neutrino emissivity in the reaction (|^ is Q ~ 10^'^ Tg 
erg cm^^ s^^. For T ~ 10^ K, the direct Urea process is about 5-6 orders of magnitude more efficient 
than the modified Urea. Therefore, sufficiently massive NSs suffer the enhanced cooling (Sect. 7). 

In a series of papers initiated by Voskresensky and Senatorov (1984, 1986) the neutrino reactions 
of the direct Urea type have been studied for the models of dense matter with highly polarized pion 
degrees of freedom. Pion condensation in this matter takes place at /> ~ po and is very efficient. Even 
for lower p, before the condensation occurs, the neutrino emissivity appears to be much higher than 
the standard one due to polarizability of pion vacuum. Cooling of NSs with this equation of state has 
been simulated recently by Schaab et al. (1997b). We will not consider these models. 

Neutrino emission due to Cooper pairing of nucleons. Onset of nucleon superfluidity 
switches on a new neutrino generation mechanism concerned with creation of Cooper pairs. The 
process had been proposed and calculated in the pioneering article by Flowers et al. (1976) and 
rediscovered independently by Voskresensky and Senatorov (1986, 1987) ten years later. Until recently. 



the process has been forgotten due to unknown reasons and has not been included into NS coohng 
simulations. It was Page (1998a) who "remembered" the process and introduced it into the cooling 
theory. Cooling simulations including this process have been performed recently by a number of 
authors, in particular, by Schaab et al. (1997b), Levenfish et al. (1998, 1999) and Yakovlev et al. (1998, 
1999). The process in question represents actually (Yakovlev et al., 1999) neutrino pair emission (any 
neutrino flavor) by a nucleon N (neutron or proton) whose dispersion relation contains an energy gap: 

N ^ N + i^ + i?. (5) 

The reaction cannot occur without superfluidity: the emission of a neutrino pair by a free nucleon 
is forbidden by energy-momentum conservation. According to Yakovlev et al. (1999) (see also Sect. 
6.1), the neutrino emissivity due to pairing of neutrons is Q^ ^ ~ lO^^Tg F(r) erg cm^^ s^^, where 
r = T/Tc, Tc is the critical temperature of superfluidity onset, and F{t) is a function, which has 
maximum F ~ 1 at r ~ 0.4. The main neutrino energy release takes place in the temperature range 
0.2 Tc < T < Tc- The efficiency of the process in this range can be compared to or even larger 
than (Sect. 6.2) the efficiency of the modified or even the direct Urea processes suppressed partly by 
the superfluidity. This determines the importance of "Cooper" neutrinos for NS cooling. Neutrino 
emission due to pairing of protons appears to be much weaker owing to the smallness of the vector 
constant of weak neutral current involving protons. 

2.3 Cooling 

Initially a NS cools mainly via the neutrino emission from its core. However, direct detection of the 
neutrino flux is possible only during a NS birth in supernova explosion. So far it was done only once, 
from the Supernova 1987A in the Large Magellanic Cloud (see, e.g., Imshennik and Nadyozhin, 1988). 
A powerful neutrino outburst in a supernova explosion lasts for about 20 s. Afterwards the neutrino 
flux decreases rapidly in time. 

Cooling of a NS is accompanied by the loss of its thermal energy which is mainly stored in the 
stellar core. The energy is carried away through two channels: first, by neutrino emission from the 
entire stellar body (mostly from the core), where the most powerful neutrino reactions take place, and 
second, by heat conduction through the internal stellar layers towards the surface, and further, by 
thermal emission of photons from the surface. The neutrino cooling dominates at the initial cooling 
stage (see below), while the photon cooling dominates later, when the neutrino luminosity becomes 
weak (it fades faster than the photon luminosity with decreasing temperature). 

From mathematical point of view, cooling simulation reduces to solving the heat diffusion equation 
within the star (Thorne, 1977) with account for the volume (neutrino emission) and surface (photon 
emission) energy sinks. As a rule, one considers one-dimensional diffusion along radial coordinate in 
a spherically-symmetric star. In about t > 10^-10^ yr, owing to the high thermal conductivity of the 
internal layers, a wide, almost isothermal region is formed within the entire core and the main fraction 
of the crust. In this case, cooling simulations are considerably simplified being reduced to solving 
the global thermal-balance equation: the loss of the thermal energy (determined by the total heat 
capacity) is governed by the neutrino and photon luminosities of the star. This approach is described 
and used in Sect. 7.2. Accordingly the main elements of the cooling theory are: (1) NS heat capacity, 
(2) neutrino luminosity, (3) dependence of the photon luminosity on internal temperature (determined 
by the thermal conductivity of the outermost stellar layers). The first two elements are discussed in 
detail in this review. 
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Figure 2: Redshifted surface temperature T° 
superfluid NS models of different mass. 



versus stellar age t according to Page and Applegate (1992) for non- 



Character of cooling depends on many parameters: equation of state of internal layers, superfluidity 
of nucleons in the stellar core, NS mass, magnetic field, chemical composition of the surface layers, 
etc. Confronting the cooling theory with observations enables one, in principle, to constrain these 
parameters. In this review we discuss the constraints on the critical temperatures of neutron and 
proton superfluidities in the NS cores. 

While simulating NS cooling one calculates the cooling curves, the dependence of the effective 
stellar surface temperature Tg on stellar age t. The effective temperature defines the photon luminosity: 
L^ = AirTi^aT^ , where a is the Stefan-Boltzmann constant. The luminosity L°°, as detected by a 
distant observer with account for the effects of General Relativity, is related to the local luminosity L^ 
as L^ = L-y {l — TZg/TZ). Here TZg = 2GM/c^ is the gravitational NS radius. The effective temperature 
of thermal radiation Tj^, detected by a distant observer, is related to the local temperature Tg by 
T- = Ts^l - TZg/lZ. 

Figure |2| shows typical cooling curves for (non-superfluid) NSs of different masses (from Page and 
Applegate, 1992, with permission of the authors). The NS models are based on the equation of state 
which opens the direct Urea process for M > 1.35 Mq. The stars with lower M undergo the standard 
cooling. 

If the direct Urea is allowed in the stellar core, the star cools much faster than that with the 
standard neutrino luminosity. Typical cooling times t ~ T/T can be estimated from thermal balance 
of a plasma element in the stellar core: t ~ CT/Q^. Here, Qy is the neutrino emissivity, and C is 
the specific heat capacity. Using the heat capacity (Sect. 3.3) of an ideal degenerate gas of neutrons 
&t p = pq and appropriate neutrino emissivity from Sect. 2.2, one obtains (see, e.g., Shapiro and 
Teukolsky, 1983; Haensel, 1994): t^^^ ~ 1 min/Tg^ , t^^^ ~ l yr /r| . It is seen that the direct Urea 
process decreases the internal stellar temperature down to 10^ K in t^ ^ ~ 1 min, and down to 10^ K 
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in several days. For the standard neutrino emissivity, the coohng time v ' is about 1 yr and 10^ yr, 
respectively. 

These estimates describe temperature variation in the stellar core at the initial cooling stage. The 
crust cools initially much slower since its own neutrino energy losses are relatively small. Sufficiently 
low thermal conductivity of the crust produces thermal insulation between the crust and the core. 
Thus, the core and the crust cool independently at the initial stage, and the NS surface radiation 
carries no information on the thermal state of the core (almost flat initial parts of the cooling curves 
in Fig. §. 

The cooling wave from the rapidly cooling core reaches the stellar surface in 10-1000 yr (see, e.g., 
Nomoto and Tsuruta, 1987; Lattimer et al., 1994), and the surface temperature Ts falls down sharply 
(Fig. ^) . In the course of the enhanced cooling Tg drops by almost one order of magnitude in several 
years, which decreases the photon luminosity by a factor of ~ 10'^. For the standard cooling, the 
temperature fall is not so significant and lasts ~ 10^ yr. 

After the thermal relaxation in the NS interiors is over the neutrino luminosity from the core 
remains higher than the photon surface luminosity. However, at this stage the core controls the surface 
temperature. This is the second relatively flat part of the cooling curves. The internal temperature 
exceeds the surface temperature by more than one order of magnitude; the main temperature gradient 
takes place in the outermost layers of the outer crust. 

Finally, a rather old star appears to be at the photon cooling stage. The neutrino luminosity 
becomes much smaller than the photon one, and the cooling curves turn out to be much steeper (Fig. 
^. The change of the cooling regime is associated with different temperature dependences of the 
neutrino and photon cooling rates. Let us remind that the neutrino emissivity is Q^ ' oc T^ for the 
direct Urea process, and Q^ "'' oc T^ for the modified Urea process (Sect. 2.2). The photon luminosity 
falls down much weaker with decreasing temperature, as ~ T^'^ (Gudmundsson et al., 1983). Simple 
estimates show that the time t^^ of transition to the photon cooling stage is mainly determined by 
the NS heat capacity and by the temperature dependence of the neutrino energy losses; typically, 
tu ~ 10^ yr. The internal temperature of the star of age t^ is mainly determined by the efficiency of 
the neutrino energy losses and by their temperature dependence. 

The appearance of superfluidity in the NS core changes drastically the cooling process. This is 
discussed in subsequent Sects. 3-8. We will study systematically the effects of superfluidity on the 
heat capacity and neutrino luminosity of NSs. This will allow us to "calibrate" theoretical cooling 
curves and use cooling NSs as "thermometers" for measuring nucleon critical temperatures in their 
cores. We will see that the nucleon superfluidity can either enhance or slow down cooling (depending 
on parameters) and it strongly reduces difference between the enhanced and standard cooling. We 
will analyse observations of thermal radiation from cooling NSs and show that the observational data 
are difficult to explain without assuming existence of superfluidity in the NS cores. 
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3 Superfluidity and heat capacity of neutron star cores 

3.1 Nucleon superfluidity 

The theory of electron superconductivity in metals was developed by Bardeen, Cooper and Schrieffer 
(1957; BCS). Superconductivity is explained by Cooper pairing of electrons caused by a weak attraction 
due to electron-phonon interaction. Superconducting state appears with decreasing temperature as 
a result of the second-order phase transition; typical critical temperatures are Tc ~ (1^10) K. For 
T < Tc, the dispersion relation of electrons contains an energy gap A, with A ~ k^Tc at T <C Tc. 

One year after the publication of the BCS theory Bohr et al. (1958) suggested that the phenomenon 
much in common with superconductivity could appear in the systems of nucleons in atomic nuclei. 
Cooper pairing of nucleons could occur due to nuclear attraction. It was expected that the gap in the 
nucleon spectrum, A ~ 1 MeV {Tc ~ 10^'^ K), was many orders of magnitude larger than for electrons 
in metals. Later the presence of pair correlations of nucleons in atomic nuclei and associated energy 
gap in the nucleon spectrum has been investigated theoretically and confirmed experimentally (Nobel 
Prize of Bohr, Mottelson and Reynoter in 1975). 

Migdal (1959) was one of the first who applied the BCS theory to atomic nuclei. He noticed also 
that neutron superfluidity caused by nuclear forces could occur in neutron matter of the inner NS 
layers where critical temperatures Tc ~ 10^*^ K could be expected. 

BCS equations which describe symmetric nuclear matter in atomic nuclei and asymmetric neutron- 
rich NS matter have much in common but have also some differences. For instance, pairing in atomic 
nuclei takes place in the singlet-state of a nucleon pair. In this case, the energy gap is isotropic, 
independent of orientation of nucleon momentum. On the other hand, one can expect triplet-state 
pairing in NS matter (see below) which leads to the anisotropic gap. 

Calculations of the energy gap in symmetric nuclear matter have been carrying out since 1959, 
starting from the classical papers by Cooper et al. (1959) and Migdal (1959). Without pretending to 
give complete description of this activity (see, e.g., Chen et al., 1993, for a review) let us mention the 
early articles by Ishihara et al. (1963), Henley and Wilets (1964), Kennedy et al. (1964), and Kennedy 
(1968) who studied the gap in the nucleon spectrum for different model potentials (Tamagaki, 1968) of 
nucleon-nucleon interaction. It was shown by the middle of 1960s that the gap was extremely sensitive 
to the repulsive part of the potential and to the effective masses of nucleons in nuclear matter. These 
conclusions are qualitatively correct for NSs (see below). 

It is well known that the BCS theory is used also to describe superfluidity in liquid ^He. The 
foundation of the theory was built in the paper by Cooper et al. (1959) mentioned above. Super- 
fluidity in ^He (caused by interatomic attraction) is quite different from that in nucleon matter; the 
critical temperature in ^He at normal pressure is as small as 2.6 mK. However, there is one important 
similarity: pairing can occur in the triplet state of interacting particles (with orbital momentum 1=1), 
which leads to anisotropic gap. The BCS equations are similar in these cases. The flrst article devoted 
to the triplet-state pairing was written by Anderson and Morel (1961). While deriving the equation 
for the anisotropic gap the authors overlooked contribution from one of three triplet spin states. This 
inaccuracy was corrected by Balian and Werthamer (1963). 

Let us outline microscopic theories of superfluidity in NSs. Five years after Migdal (1959), Ginzburg 
and Kirzhnits (1964) published a brief article where they estimated the gap produced by the singlet- 
state pairing of neutrons at density p = lO^'^^-lO^^ g cm~'^ and obtained A ~(5-20) MeV. Very serious 
step was made by Wolf (1966). He showed that the singlet-state neutron pairing takes place in the 
inner NS crust (p < po), but disappears in the core, since the singlet-state nn-interaction becomes 
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repulsive in high-density matter. The number density of protons in the NS core is much smaller than 
the number density of neutrons (Sect. 2.1). Accordingly the single-state pp-interaction is attractive 
there and leads to the proton pairing. 

Before the discovery of pulsars in 1967 the theory of superfluidity in NSs was developing rather 
slowly. A detailed review of "pre-pulsar" articles was given by Ginzburg (1969). The discovery of 
pulsars initiated a great interest to the theory. Baym et al. (1969) analyzed macroscopic consequences 
of the neutron superfluidity (rotation of superfluid component in the form of quantized vortices) and 
proton superfluidity - superconductivity (splitting of the internal stellar magnetic field into Abrikosov 
magnetic flux tubes). This paper made foundation of the modern theory which explained pulsar 
glitches by interaction of normal and superfluid components of matter inside NSs (see, e.g., Pines, 
1991). 

Very important contribution to the theory was made by Hoffberg et al. (1970). They showed 
that the triplet-state ^P2 interaction of neutrons at p > po was attractive. Thus, the triplet-state 
neutron superfluidity with anisotropic gap can occur in the NS core. The authors performed the first 
calculations of the critical temperature of the triplet-state neutron superfluidity in NSs. 

The article by Hoffberg et al. (1970) was followed by many others, where the theory was developed 
and nucleon critical temperatures in NSs were calculated. Superfluidities of various types have been 
considered using different model potentials of nucleon-nucleon interaction. Some authors have made 
use of the potentials renormalized with account for polarization properties of matter on the basis 
of various many-body theories. A comparative analysis of different approaches has been done, for 
instance, by Amundsen and 0stgaard (1985a, b), Wambach et al. (1991), Baldo et al. (1992), and 
Chen et al. (1993). A detailed review has been written by Takatsuka and Tamagaki (1993) not long 
ago. Many articles have been devoted to the triplet-state superfluidity of neutrons in the NS cores 
(see, e.g., Hoffberg et al., 1970; Tamagaki, 1970; Takatsuka and Tamagaki, 1971, 1993, 1997a, b, c; 
Takatsuka, 1972a, 1972b; Amundsen and 0stgaard, 1985b; Baldo et al., 1992; Elgar0y et al., 1996a, 
c). Let us mention model calculations by Muzikar et al. (1980) of the triplet-state neutron pairing 
in magnetized NS cores. According to these calculations, the magnetic field B > 10^^ G makes 
the triplet-state superfluidity with the nodes at the Fermi-surface energetically more preferable than 
the familiar triplet-state superfluidity without nodes (Sect. 3.2). Singlet-state proton superfluidity in 
the NS cores has been considered thoroughly by many authors (e.g., Chao et al., 1972; Takatsuka, 
1973; Amundsen and 0stgaard, 1985a; Baldo et al., 1990, 1992; Chen et al., 1993; Takatsuka and 
Tamagaki, 1993, 1997c; Elgar0y et al., 1996a, b). Elgar0y and Hjorth-Jensen (1998) have studied 
a singlet-state proton superfluidity in the NS cores and estimated maximum densities at which this 
superfluidity disappears. Let us mention also the articles by Sedrakian et al. (1997) and Civitareze 
et al. (1997), where the possibility of neutron-proton Cooper pairing in uniform nucleon matter has 
been analyzed. This pairing is possible in symmetric nuclear matter but does not occur in NSs due 
to large difference of neutron and proton Fermi momenta. Many authors (e.g., Hoffberg et al., 1970; 
Krotscheck, 1972; Clark et al., 1976; Takatsuka, 1984; Amundsen and 0stgaard, 1985a; Chen et al., 
1986, 1993; Ainsworth et al., 1989; Baldo et al., 1990, 1992; Wambach et al., 1991, 1993; Takatsuka 
and Tamagaki, 1993; Brogha et al., 1994; Elgar0y et al., 1996b, 1996d; Pethick and Ravenhall, 1998) 
have considered the singlet-state pairing of free neutrons in the inner NS crust. Nucleons within atomic 
nuclei in the inner crust can also suffer pairing which is, however, weaker than for free neutrons (see, 
e.g., Takatsuka, 1984; De Blasio and Lazzari, 1995; Elgar0y et al., 1996d; Barranco et al., 1998; Pethick 
and Ravenhall, 1998). 

Although we will not consider exotic models of NS cores (Sect. 2.1), let us mention that superfluidity 
is possible in these models as well. For instance, Takatsuka and Tamagaki (1993) have reviewed 
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calculations of neutron and proton superfluid gaps in pion condensed matter. They have obtained 
also new results in this field (Takatsuka and Tamagaki, 1997a, b, c). The same authors (Takatsuka 
and Tamagaki, 1995) have studied the nucleon superfluidity in the presence of kaon condensation. 
The pion or kaon condensates mix strongly neutron and proton states which may induce triplet-state 
pairing of quasi-protons. Some authors have discussed superfluidity in quark matter (e.g., Balian and 
Love, 1984; Iwasaki, 1995; Schaab et al., 1996, 1997a). If hyperons appear in npe-matter, they can 
also be in a superfluid state (Balberg and Barnea, 1998). Electrons and muons, which interact via 
Coulomb forces, can be, in principle, superfluid as well. However, corresponding critical temperatures 
are too low to be of practical interest. For instance, according to an estimate by Ginzburg (1969) the 
critical temperature of degenerate electrons at p > 100 g cm~'^ does not exceed 1 K, while the internal 
temperature of a cooling NS of age t < 10^ does not fall below 10® (Sect. 7). 

Let us summarize properties of nucleon superfluidity in the cores and crusts of NSs with the 
standard {npe) composition. 

(1) Singlet-state neutron superfluidity exists in the inner NS crust and disappears at the density p ~ 
po, at which an effective neutron-neutron singlet-state attraction transforms into repulsion. Density 
dependence of the critical temperature Ten has maximum at a subnuclear density. Maximum values 
of Ten range from 10® K to 10^^ K, for different models of dense matter. 

(2) Triplet-state neutron superfluidity appears in the NS core at p > po owing to an effective 
attraction between neutrons in the triplet state. Density dependence of the appropriate critical tem- 
perature, as a rule, has maximum at supranuclear density. Maximum values of Ten vary from 10*^ K 
to 10^" K for different microscopic models. 

(3) Protons in the NS cores can suffer singlet-state pairing. The dependence Tcp{p) has usually 
maximum at a supranuclear density, and the maximum values of Tcp range from 10*^ K to lO^'' K for 
different models of matter. 

(4) The critical temperature is very sensitive to the strength of repulsive core of nucleon-nucleon 
interaction (see, e.g., Takatsuka and Tamagaki, 1971, 1993; Krotscheck, 1972). Tc increases strongly 
with core softening, i.e., with increasing attraction between nucleons. Even a weak additional attrac- 
tion (e.g., due to inclusion of coupling between the ^P2 and ^F2 states for the triplet-state pairing of 
neutrons: see Takatsuka, 1972a; Takatsuka and Tamagaki, 1993) may increase Tc by several orders of 
magnitude. Generally, the superfluidity is stronger for softer equation of state. 

(5) The critical temperature falls down rapidly with decreasing the effective mass of nucleons in 
dense matter (see, e.g., Takatsuka and Tamagaki, 1971, 1993; Chao et al., 1993), i.e., with decreasing 
density of state of nucleons at the Fermi surface. If the effective mass is sufficiently low (m^ < 0.5mj\[) 
the superfluidity may entirely disappear. 

(6) The critical temperature depends strongly on the method of description of many-body (in- 
medium) effects. For instance, the in-medium effects for neutrons in the NS crusts decrease Ten typi- 
cally by several times (e.g., Clark et al., 1976; Wambach et al., 1991, 1993; Takatsuka and Tamagaki, 
1993). 

For illustration, in Fig. we present the critical temperature of the triplet-state neutron superflu- 
idity and the singlet-state proton superfluidity obtained by different authors. 

The HGRR curve has been calculated by Hoffberg et al. (1970) for neutron superfluidity using the 
Tabakin (1964) model of nucleon-nucleon potential which reproduces quite well experimental phases of 
nucleon scattering at energies < 320 keV. Hoffberg et al. (1970) have neglected the in- medium effects; 
the effective neutron mass has been set equal to its bare mass. The TT-curve has been calculated 
by Takatsuka and Tamagaki (1971) for the neutron superfluidity under the same assumptions but 
using the one-pion-exchange model of nucleon interaction (OPEC ^0-1) with somewhat harder core. 
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Figure 3: Critical temperatures of triplet-state neutron superfluidity (solid lines) and singlet-state proton superfluidity 
(dashed lined) versus density according to different authors: HGRR — Hoffberg et al. (1970), TT — Takatsuka and 
Tamagaki (1971), CCY — Chao et al. (1972), T — Takatsuka (1973), AO — Amundsen and 0stgaard (1985a, b), BCLL 
— Baldo et al. (1992) (see text for details). 

This lowers Ten- The sohd AO-curve has been obtained by Amundsen and 0stgaard (1985b) for the 
neutron superfluidity using similar one-pion-exchange approach (OPEG), but the effective neutron 
mass has been determined self-consistently and has appeared to be lower than the bare neutron mass. 
This additionally lowers Ten- The BCLL-curve is a result of recent calculations by Baldo et al. (1992) 
using an empirical nucleon potential (Argonne ^14) which reproduces accurately laboratory data on 
nucleon scattering. The potential and neutron mass have not been renormalized to account for the 
in-medium effects. As a result, the values of Ten are high again. 

The dashed CCY curve shows the results by Chao et al. (1972) for proton superfluidity using 
the simplified Serber potential which describes satisfactorily the experimental phases of singlet-state 
nucleon-nucleon scattering. The effective proton mass has been determined self-consistently. The T 
curve presents the results by Takatsuka (1973) obtained under the same assumption but using a more 
realistic OPEG potential. The dashed AO curve has been calculated by Amundsen and Ostgaard 
(1985a) using the same one-pion-exchange approach but another method of evaluation of the effective 
proton mass which leads to noticeably different critical temperatures Tcp. Finally, the dashed BCLL 
curve has been derived by Baldo et al. (1992) under the same assumptions as the solid BCLL curve. 
A neglect of polarization effects leads to the high critical temperature Tcp. 

Let us emphasize large scatter of Ten and Tcp. For instance, the BCLL curves (Baldo et al., 
1992) mentioned above give very high critical temperatures in matter of supranuclear density while, 
for instance, calculations by Takatsuka and Tamagaki (1997b) indicate significant decrease of critical 
temperatures at densities of several po, at which the direct Urea process is open. Under these condi- 
tions, while analysing the effects of superfluidity, we will not use any speciflc microscopic results but 
will treat Ten and Tcp as free parameters. 
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Table 1: Three superfluidity types of study 





Superfluidity type 


A 


F{d) 


kBTjA{0) 


A 


^So 


1 


1 


0.5669 


B 


3P2 {mj = 0) 


1/2 


(l + Scos^t?) 


0.8416 


C 


3p2(|mj| = 2) 


3/2 


sin •& 


0.4926 



3.2 Energy gaps and critical temperatures 

Thus, the cases of the ^Sq or ^P2 pairing in the NS cores are of special interest. The ^P2 pairing 
in npe matter occurs mainly in the system of neutrons. While studying this pairing one should take 
into account the states with different projections mj of the total pair momentum onto quantization 
axis: \mj\ = 0, 1, 2. The actual (energetically favourable) state may be a superposition of states with 
different mj (see, e.g., Amundsen and 0stgaard, 1985b; Baldo et al., 1992). Owing to uncertainties 
of microscopic theories this state is still unknown; it depends possibly on density and temperature. In 
simulations of NS cooling, one usually considers the triplet-state pairing with mj = (excluding the 
recent article by Schaab et al., 1998b). Below we will consider the ^P2 superfluidities with ttij = 
and \mj\ = 2, since their effects on the heat capacity and NS neutrino luminosity are qualitatively 
different. 

For certainty, we will analyze the BCS superfluidity for an ensemble of almost free (quasi-)nucle- 
ons. The superfluidity types ^Sq, ^P2 ("^j = 0) and ^P2 {rrij = 2) under study will be denoted by A, 
B and C, respectively (Table ||). Superfluidity onset is accompanied by the appearance of the energy 
gap 5 in momentum dependence of the nucleon energy e(p). Near the Fermi surface (|p — Pp\ *C Pp) 
one has (e.g., Lifshitz and Pitaevskii, 1980) 



M 



Y 5^ + T]"^ at p < pp] e = fi+ \l 5^ + rf' at p> pp. 



(6) 



Here, tj = 'Vpip — Pp), Vp and pp are the nucleon Fermi velocity and momentum, respectively, n is the 
chemical potential; it is assumed that 6 <^ n. In the cases of study, S'^ = A^(T)F (■!?), where A(r) 
is an amplitude, which determines temperature dependence of the gap, F^i!}) describes dependence of 
6 on angle "& between the quantization axis and particle momentum p. The quantities A and F are 
determined by the superfluidity type (Table |lj). In case A the gap is isotropic, 6 = A. In cases B 
and C the gap is anisotropic (depends on i?). Let us notice that in case C the gap vanishes at the 
Fermi-sphere poles at any temperature (since F{0) = F{tt) = 0), i.e., the superfluidity does not affect 
the nucleons which move along the quantization axis. 

The gap amplitude A(r) is determined by the BCS equation (see, e.g., Lifshitz and Pitaevskii, 
1980; Tamagaki 1970) which can be written as 



In 
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47r Jo z 



(7) 



where Aq = A(0), dfi is a solid angle element in the direction of p, / = (1 + e^) ^ is the Fermi-Dirac 
distribution, A is a numeric coefficient (Table |l]). 
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Using Eq. (^) one can easily obtain the values of k-QTc/ ^o presented in Table ||. It is convenient to 

introduce the variables 

A(r) T 

kbT Tc 

The dimensionless gap amplitude v describes the temperature dependence of the gap. It is determined 
by the superfluidity type and the dimensionless temperature r. In case A the amplitude v corresponds 
to the isotropic gap, in case B it corresponds to the minimum and in case C to the maximum gap at 
the Fermi surface. In these notations, the dimensionless gap y has the form: 

yA=VA, yB = % Vl + Scos^^, yc = Vc sin 'd. (10) 

Using Eq. (^) one can obtain the asymptotes of the gap amplitude near the critical temperature and 
in the limit of the so called strong superfluidity (T <^ Tc). For instance, at T — > Tc, T < Tc (r — > 1) 
we have (see, e.g., Lifshitz and Pitaevskii, 1980; Levenfish and Yakovlev, 1994a): v = /3\/l — r, where 
Pa = 3.063, f3B = 1-977, Pc = 3.425. For T <C T^ we obtain v = Ao/{kBTcT). Levenfish and Yakovlev 
(1993, 1994a) calculated v = v{t) for intermediate r and proposed the analytic fits of the numeric 
data: 

r. f. .rn 0-157 1.764 

V, = VI -r 1.456 ^ H 

V ^Jt T 

vj^ = Vr^T (^0.7893 + h}^\ , (11) 



Vc = — (2.030 - 0.4903t^ + 0.1727r^ 

r ^ 



These fits will be useful for evaluating the heat capacity and neutrino luminosity in the superfluid NS 
cores. 

Analytic fits presented here and below reproduce numeric data with mean errors about 1-2%, while 
the maximum error does not exceed 5%. This fit accuracy is more than sufficient for implications to 
NS cooling. Our fits reproduce also various asymptotes of the quantities of study. 

Notice that for the triplet-state superfluidity with \mj\ = 1 or for a triplet-state superfluidity 
described by superposition of states with different rrij, the anisotropic gap 5 depends not only on iD, 
but also on azimuthal angle (p of nucleon momentum p. A study of the effects of such superfluidity 
on the heat capacity and neutrino luminosity is complicated and has not been performed yet. 

3.3 Heat capacity of superfluid neutron star cores 

Consider npe matter of the NS cores. The specific (per unit volume) heat capacity is equal to the sum 
of partial heat capacities of all particle species j: C = J2i Cj, where j = n, p, e. Since the particles are 
strongly degenerate, the heat capacities at constant volume and pressure are nearly the same (Lifshitz 
and Pitaevskii, 1980) and we make no difference between them. A partial heat capacity is determined 
by the standard thermodynamic formula (Lifshitz and Pitaevskii, 1980) 

where £j and pj stand for the particle energy and momentum, respectively, /Xj is the chemical potential, 
and fj is the Fermi-Dirac distribution. 
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In our case, the electrons constitute an almost ideal, strongly degenerate, ultra-relativistic gas. 

Accordingly, 

^ mlVf^^kVr iQ/"-e\^^^^ erg , , 

C^e = ^^3^-5.67x1019 — Tg ^-, 13 

3?i'^ \noJ cm'^K 

where m* = ^e/c^ ~ Ppe/*^' PFe ^^*i "-e denote, respectively, the electron Fermi momentum and number 
density, no = 0.16 fm~^ is the standard nucleon number density in atomic nuclei. In the presence of 
muons, one should add their partial heat capacity C^ to the total heat capacity; C^ is similar to Ce 
with the only difference that muons may be non-relativistic. 

Neutrons and protons {j = N = n and p) in the NS cores constitute a non-relativistic, strongly 
non-ideal Fermi-liquid (Sect. 2.1). The partial heat capacity (12) of normal nucleons is: 



^^H^ - 1-61 X 10^° ^ f^V^' T, -^ . 
3?i'^ rriN Kuq J cm'^ K 



C^^ = -Nt^m^'^B- ^ l_gl ^ 1020 !^ ( !1^ ) Tg ^^^. (14) 



Here, pp^ is the nucleon Fermi momentum determined by the nucleon number density n^v, and m^ is 
the effective mass of the (quasi-)nucleon in dense matter. Notice that the main contribution into the 
heat capacity comes from nucleons with energies near the Fermi level, \£n — /^Af| ^ k-^T, which can 
participate in processes of energy exchange ~ k-oT. 

For T < Tc, Eq. (|l^) should include the energy gap (P). Generally, the nucleon heat capacity can 
be written as 

CN = CmR{T), (15) 



where Cnq is the partial heat capacity of normal nucleons (14), and the factor R describes variation 
of the heat capacity by the superfluidity. R depends on superfluidity type and on dimensionless 
temperature r. Clearly, i? (T) = 1 for T > Tc (r > 1). 

A general expression for R is obtained from Eq. ( [l2|) by introducing the dispersion relation (gj), 
dimensionless variables (S) and by taking into account Eq. (p!^): 



Notice that the Fermi-Dirac distribution / depends on T directly as well as functionally, through the 
dispersion relation (|6|). Therefore, d//dT contains derivative of the gap amplitude A(T) with respect 
toT. 

The heat capacity of superfluid Fermi-systems is well known in the physical literature (especially 
for the singlet-state pairing) but not in the astrophysics of NSs. The heat capacity Ca(T) for case 
A has been calculated by Miihlschlegel (1959) and is described in textbooks (Lifshitz and Pitaevskii, 
1980). Maxwell (1979) proposed a fit of C^iT) in the temperature range 0.2 Tc <T< Tc- Anderson 
and Morel (1961) derived an asymptote of CciT) at T <C Tc for superfiuidity of type C. Simple 
expressions for i?^, Rb and Re convenient for evaluating the nucleon heat capacity in the NS cores 
were obtained by Levenfish and Yakovlev (1994a) and are given below. 

When temperature falls below the critical value, the heat capacity suffers a jump produced by 
latent heat release at the phase transition. In case A the jump is given by Ra{Tc) = 2.426 (Lifshitz 
and Pitaevskii, 1980), while in cases B and C the jump is Rb{Tc) = Rc{Tc) = 2.188. These jumps 
affect NS cooling (Sect. 7.2). 

The energy gap in the nucleon spectrum strongly reduces the heat capacity for temperatures much 
below the critical one (T <C Tc, v — > oo). Let us remind, that the main contribution into the heat 
capacity in non-superfluid matter would come from nucleons with energies |e — ^| < k-oT. However, 
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Figure 4: Factors 7?, which describe variation of the heat capacity by superfluidity, versus T/Tc- Letters near the 
curves show superfluidity type. Dashed curves A* and B* present simplified factors K* for superfluidities A and B, 
which were commonly used in applications earlier (see text). 

in the low-temperature limit T <^ Tc for all nucleons in cases A or B and almost for all nucleons 
(excluding a small fraction of particles near the poles of the Fermi sphere) in case C the gap S is much 
larger than k^T. Therefore, the "heat-capacious" nucleons with energies |e — /i| < fceT are either 
absent or almost absent, which suppresses the heat capacity. Accordingly, we have the asymptotes: 



i?^ = X2^5/2^_. 



3.149 



7r3/2 
vr 



Rr, = ^^v^ e"^ 



r5/2 
0.7781 



exp 



T^ 



exp 



1.764 



(17) 



Re 



111 
5v^ 



3.353 r^. 



Here and hereafter we set f = f^, f = f^ or t; = vc in the factors R, labeled by indices A, B or C, 
respectively. With decreasing temperature, Ra and Rb are reduced exponentially, while Re is reduced 
much slower, proportionally to T^. The exponential reduction is associated with exponentially small 
"efficiency of excitation" of quasi- nucleons near the Fermi surface in the presence of a large gap. In 
the case of the heat capacity, this effect strongly reduces "heat exchange", but, generally, all nucleon 
reactions (transitions) are reduced. A small exponent is absent in case C because the gap 5c = 5c {T, i?) 
vanishes for nucleons at the Fermi-sphere poles. The superfluidity disappears near the poles, and these 
nucleons have "normal" heat capacity. However the fraction of these nucleons decreases proportionally 
to T^, which leads to a power-law reduction of the heat capacity of the nucleon liquid as a whole. 
The results of numeric calculations of R, as well as the asymptotes for v ^ and v ^ oo are 
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Figure 5: Factors 7?4 (curve A) and Rj^ (curve B), as in Fig. W, but versus dimensionless gap amplitude v. 
described by the simple fit expressions (Levenfish and Yakovlev, 1994a) 

Rji = ( 0.4186 + y'(1.007)2 + (0.5010 t;)2j exp M .456 - ^^(1.456)2 + i;2 J , 
Rb = ( 0.6893 + ^^(0.790)2 + (0.2824 7;)2j exp M .934 - y^(1.934)2 + z;2 
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Equations (p^), combined with (|ll|), enable one to evaluate easily the factor i? as a function of the 
dimensionless temperature r. 

The dependence of i? on r for superfluidity types A, B, C is shown in Fig. ^. The factors Ra and 
Rb as functions of the dimensionless gap amplitude v are close which is illustrated in Fig. |5l 

In simulations of the NS cooling before publication of the articles by Levenfish and Yakovlev (1993, 
1994a), as far as we know, one considered the superfluidity of A or B type and used the simplifled 
factors of the form R* = exp(-i;*), where v\ = dA/ikjiT) = 1.764/r, v*b = 5max/(feBT) = 2.376/r, 
^max = 2Ab(0) is the maximum energy gap (5b(0, ??) as a function of ?? at T = 0. These factors were 
introduced from physical consideration (see, e.g., Page and Applegate, 1992) basing on the suggestion 
by Maxwell (1979). As seen from Fig. ^ the simplified factors do not reproduce the heat capacity jump 
at T = Tc and strongly overestimate reduction of the heat capacity at T <C Tc. Indeed, the correct 
asymptotes in the limit T <^ Tc contain large pre-exponents omitted in the simplified expressions. 
In addition, the factors Rb and Rb differ by exponent arguments. Instead of the correct argument 
vb = Smm/{kBT) (see Eq. (|T7|)) the simplified factor contains the argument Vb = 5max/(^B^) = 
2^111111/(^6^), which is twice larger. As a result, for instance, at T = 0.1 Tc the exact factor Ra exceeds 
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the simplified factor i?^ by about three orders of magnitude, while the exact factor Rb exceeds -Rg 
by seven orders of magnitude. 

Although the local heat capacities of particles are determined by a given equation of state, the 
relative contribution of n, p and e into the total heat capacity of a non-superfluid NS core does not 
depend strongly on the equation of state (Page, 1994). The neutron heat capacity is about 3/4 of 
the total one, the proton heat capacity is about 1/4 of the total one, and the electron contribution is 
only about 5%. Therefore, a strong superfluidity of n (with normal p) reduces the NS heat capacity 
by a factor of ~ 4. In the opposite case (strong proton superfluidity and normal neutrons) the total 
heat capacity is reduced only by ~ 25%. If neutrons and protons are strongly superfluid at once, only 
~ 5% of the heat capacity of the normal core survives. In Sect. 7.2 we will see that the asymmetry of 
the nucleon composition of the NS cores affects the NS cooling. 
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4 Reduction of direct Urea process by nucleon superfluidity 



4.1 Direct Urea proeess in non-superfluid matter 

Neutrino cooling of a NS with the standard nuclear composition is determined by the direct Urea 
process (B) involving neutrons and protons; the process is allowed (Lattimer et al., 1991) for many 
realistic equations of state in the inner NS core. In the outer core, the process is forbidden by 
momentum conservation of reacting particles (Sect. 2.2). 

The process (Q) consists of the direct and inverse reactions. Under beta-equilibrium condition, their 
rates are equal, and it is sufficient to calculate the neutrino emissivity in the first reaction and double 
the result. The total emissivity will be labeled by the upperscript (D); in a normal (non-superfluid) 
matter, according to Lattimer et al. (1991), it is given by (h = c=kB = l)'- 
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where pj is a nucleon momentum (j=l corresponds to neutron and j=2 to proton), Pe and Eg are 
electron momentum and energy, respectively, pjy and e^ are antineutrino momentum and energy, the 
delta function S{Ef — Ei) describes energy conservation, and (^(P/ — Pj) describes conservation of 
momentum of particles in the initial (f) and final (/) states. Furthermore, C stands for the product 
of Fermi-Dirac distributions or corresponding blocking factors of nucleons and electron, |Mp is the 
squared reaction amplitude. Summation is over initial and final spin states. The step function Q 
forbids the reaction in matter of not too high density (Sect. 2.2): B = 1 if J'f„, p-pp^ PPe satisfy the 
"triangle condition" and = otherwise. 

Nucleons and electrons in the NS core are strongly degenerate and the main contribution into the 
integral ( [l9| ) comes from narrow regions of momentum space near the corresponding Fermi surfaces. 
Thus one can set p = pp in all smooth functions under the integral. For non-relativistic nucleons, 
the quantity |Mp is independent of particle momenta and can be taken out of the integral. The 
remaining integral is evaluated by the standard method of decomposition of integration over directions 
and magnitudes of particle momenta. Introducing the dimensionless variables in accordance with Eq. 
(|8|), the neutrino emissivity (p!9|) can be written as (e.g., Levenfish and Yakovlev, 1994b) 
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The quantity A contains integrals over orientations of particle momenta {j = 1, 2, 3, and 4 corresponds 
to n, p, e, and v, respectively); dOj is solid angle element in direction of pj. All vector lengths pj 
in the delta-function must be set equal to the corresponding Fermi momenta. The quantity /, given 
by Eq. (P2|) , includes integration over dimensionless energies of neutrino x^, = p^/T = £y/T and other 
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particles Xj = VpAp — pp-)/T (see Eq. (^); fj = [exp(xj) + 1]~^ is the Fermi-Dirac distribution. For 
particles j=2 and 3, we have transformed [1 — /(x)] -^ fix) by replacing x — > —x. Finally, S contains 
products of the density states at the Fermi surfaces of n, p and e; m* is an effective particle mass, 
with m* = ^e/c^- As shown, for instance, by Prakash et al. (1992) the quantity |Mp, summed over 
spins of final particles and averaged over spin of initial neutron and over directions of electron and 
neutrino momenta is Gp cos^^c {fy + ^g'^), where Gp = 1.436 x 10~^^ erg cm^ is the Fermi constant 
of weak interaction; /y = 1 and qa = 1-26 are, respectively, the vector and axial-vector constants for 
the reaction of study, and 9c is the Cabibbo angle (sin^c = 0.231). 

In the absence of superfluidity, the integrals A and / are standard (e.g., Shapiro and Teukolsky, 
1983): ^0 = 327r^/ (pp^pp »pg) and Iq = 4577r^/5040. Thus, the neutrino emissivity of the direct 
Urea process in non-superfluid matter (in the standard physical units) is (e.g., Prakash et al., 1992) 

iTi^ 457 vr o o o n m*m*m* ^ 

(D) ^ ^O^ 2 2 Q (^2 ^ 3 2) n p e ^^^^6 q 

^° 10080 ^ ^^ ^' n^^(? 

/« \V3 rn* m* ^ ' 

«. 4.00 X lo^M ^ ) ^^r|e ^^. 

Vno/ m^ cm'^-s 

Here, as before, no = 0.16 fm~^. 

Similar formula is valid for the direct Urea processes involving hyperons. Notice that /y and 
QA can be renormalized under the action of in-medium effects. Exact calculation of these effects is 
complicated. Here and below, for certainty, we will use non-renormalized values. Notice also that the 
direct Urea process is affected by the strong magnetic field, B > 10^^ G (Baiko and Yakovlev, 1999). 

4.2 Direct Urea proeess in superfluid matter 

The main contribution into neutrino generation in non-superfluid matter comes from nucleons with 
energies |e — ^| < k^T. The nucleon spectrum in superfluid matter contains the energy gap which 
suppresses the reaction. The essence of suppression is the same as in the case of suppression of heat 
capacity by strong superfluidity (Sect. 3.3). In contrast to the heat capacity, even a weak superfluidity 
usually suppresses the neutrino reactions (excluding the neutrino emission due to Cooper pairing of 
nucleons. Sects. 2.2, 6.1). 

The expression for the neutrino emissivity in the direct Urea process can be generalized to the 
case of superfluid matter by introducing the energy gap into the e(pj) dependence (see Eq. (^)) for 
nucleons. Let protons (j = 2) suffer Cooper pairing of type A, while neutrons (j = 1) suffer pairing 
of types A, B or C. In order to incorporate superfluidity in Eq. (^) it is sufficient to replace Xj — > Zj 
(see Eq. (^)) and introduce averaging over orientations of pi. Then the emissivity can be written as 

Q(D)^g(D)^(D)^ i?(°)(^l,«2) = Y^=j^J^y^^y^) =/o"^'^^ Sin^ J(yi,y2), (25) 

where y2 = V2, yi is given by Eq. ([To|), 

1 ("+00 ^ r+oo r+oo /■+oo 

J{yi,y2) = Y / dXiy^l / dXi/(zi) / dX2 7(2:2) / dXef{Xe)d{Xu-Zi-Z2-Xe), (26) 

^0 Jo J—OD J— 00 J —00 

dQ is solid angle element along pi, t? is an angle between pi and the quantization axis; Qq is the 
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emissivity (|24| ) in non-superfluid matter; and R^ ' is the factor which describes superfluid suppression 
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of the reaction. One formally has R^ ' = 1 for normal nucleons, while i?^ ' < 1 in the presence 
of superfluidity. The integrals ([25| ) and (^) have been calculated by Levenfish and Yakovlev (1993, 
1994b) for different combinations of neutron and proton superfluidities. The results are discussed 
below. 

4.3 Superfluidity of neutrons or protons 

Consideration of the cases of superfluid neutrons or protons is similar. For instance, let neutrons be 
superfluid. Then we can set Z2 = X2 in Eqs. ( p5[ ) and (p6|), and w ' will depend on one argument 
vi = V and on superfluidity type. For r = T/Tc > 1, as mentioned above, w ' = 1. For a strong 
superfluidity (r <C 1, f >» 1) the neutrino emission is greatly suppressed {R^ ^ ^ 1). The asymptotes 
of R^ ' in the cases of strong superfluidity of types A, B and C read: 



„(D) 252 /^ 5 5 , , 0.0163 / 1.764\ 

^^ 457^ V 2 ^ ■ '^P^-^^ = -^^ '^P [~~r) ' 

„(D) 126 ^ . ^ 0.00123 / 1.188\ 

Rr = r^ ■^ exp -V = ^- exp , 27 

^ 4577r5^/3 ^ ^ t5 "^ \ t J ' ^ ^ 

(D) 6029 vr^ „ 

^c = ..o. 2 = 2.634 r^. 
^ 5484 f 2 

The asymptotes of Rj^ and -Rg are seen to be exponential while the asymptote of R^ is power-law. 
The latter circumstance is associated with the nodes of the gap at t? = and vr. Comparing with the 
results of Sect. 3.3, we see that these asymptotes are similar to the asymptotes of the nucleon heat 
capacity in the limit of strong superfluidity. 

The asymptotes (p^) and numeric values of the integrals (^) and (2£) for intermediate v can be 
fitted by: 
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exp 2.701 - V(2.701)2 + v^ ] , (28) 



^(B) 0.5 + (0.09226^)2 1 / / 

Rr = } ^:t5 ; TT + - exp 1 - Jl + 0.4129 u 

^ l + (0.1821z;)2 + (0.16736^)4 2 ^V ^ 

Equations (p8|), with account for (0), enable one to calculate R^ ' for any r. The results are presented 
in Fig. |6[ 

In simulations of NS cooling prior to publication of the above results (Levenfish and Yakovlev, 1993, 
1994b) one usually used (e.g.. Page and Applegate 1992) the simplified factors of reduction of the direct 
Urea process by the superfluidity of type A or B. These factors were proposed by Maxwell (1979) 
from simple (and inadequate) consideration. Comparison of the accurate and simplified factors shows 
that the latter factors strongly overestimate the direct Urea reduction. For instance, for T = 0.1 Tc 
the accurate factor R)^ appears to be about four orders of magnitude larger than the simplified one, 
and the accurate factor R^ is more than seven orders of magnitude larger. The difference has the 
same nature as the difference of accurate and simplified heat-capacity factors discussed in Sect. 3.3. 
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Figure 6: Superfluid reduction factors of the direct Urea process versus r = T/T^. Letters near the curves indicate 
superfluidity type (see Table 1): A, B and C mean superfluidity of neutrons or protons, while AA, BA and CA mean 
superfluidity of neutrons and protons. 

4.4 Superfluidity of neutrons and protons 

Let us consider reduction of the direct Urea process by combined actions of superfluidities of neutrons 
and protons of types (AA), (BA) and (CA). In these cases, the reduction factor R^ ' depends on 
two arguments, vi = Vn and V2 = Vp. 

Calculation of w ' in the presence of the neutron and proton superfluidities is complicated; let 
us consider the factor R)^^, as an example. According to Eqs. (^) and (25), for the singlet-state 
pairing R\J{vi,V2) = J{vi,V2) = J(f2,fi), where yi = vi, y2 = V2- Clearly, R\J{0,0) = 1. If both 
superfluidities are strong {vi » 1 and V2 ^ 1) and f2 — ui » \/v2, the asymptote of the reduction 
factor is: 
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where s = 
asymptote 



— vf. In the limit vi <C V2 Eq. (|30| ) gives K = vf/^O' which corresponds to the 
271) of Rj^ . In another limit 



V2 ^1^2 — ^1 ^ ^2 we obtain K = (2/315) s/^l- The 
< ^/v2- One can show (Levenfish and Yakovlev, 1994b) 



asymptote (30) fails in the range | V2 
that K ~ y/v2 for vi = f 2 . 

Let us present also the fit expression which reproduces the asymptote ( p9|) and numeric values of 
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2p = M + 12.421 + Vw2 + 16.350u + 45.171, 
2g = n + 12.421 - Vw2 + 16.350^ + 45.171, 
2p, = u + v/u^2_^ 5524.8 n + 6.7737, 



2pe = u + 0.43847 + \/u;2 + 8.3680 w + 491.32, 
D = 1.52 (uiU2)^/^ (ui + U2) exp(-ui - ti2). 



m = 1.8091 + Jvf + (2.2476)2, 



n2 = 1.8091 + y^i;| + (2.2476)2, (32) 

with u = vf + V2 and w = V2 — vf. For ^2 = 0, the factor R^J^{vi, 0) agrees with the factor R}^ (vi) 
given by Eq. (p8|). 

Figure shows lines of i?^y( = const versus ri = T/Td and r2 = T/Tc2- The lines exhibit superfluid 
reduction of the direct Urea for any T, Td and Tc2. The behaviour of R\/^ at r^ + r| ^ 1 (both 
superfluidities are strong) is of special interest. In this case one can obtain an approximate relationship 
(Lattimer et al., 1991; Levenfish and Yakovlev, 1994b) 

ii(? ~ min (r?\ Rf^) , (33) 



where Ri and i?2 ^-^^ the reduction factors for superfluidities of one type. According to (33), the 
factor R\2 is mainly determined by a stronger superfluidity (1 or 2). The presence of the second, 
weaker superfluidity reduces naturally R12 , but not to a great extent. This is confirmed by the 
asymptote (|29| ) . A transition from the regime in which R\2 ~ -Ri to the regime in which R\2 ~ i?2 
takes place in the region V2 ~ i^i, where the asymptote -R12 requires special consideration. 

Using Eq. (|25|), it is not difficult to evaluate R^ ' for cases (BA) and (CA), in which the protons 
suffer singlet-state pairing while the neutrons suffer triplet-state pairing. The evaluation is reduced 
to one-dimensional integration in Eq. (p^) of the quantity J(yi, 2/2)5 fitted by Eqs. ( |3l| ) and (32). The 
results of calculations of R^j^ and R(jj^ for any T, Td and Tc2 are shown in Fig. 0. The dependence 
of R^j^ and R(j]^ on ri and T2 is similar to the dependence of Rj^j^, shown in the same figure, but 
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now R^ (Ti,r2) 7^ R^ (■^2)''"i)- The approximate expression ( p3| ) in these cases is valid as weh 
as confirmed by the asymptotes (Levenfish and Yakovlev, 1994b). Since the superfluidity of type 
C reduces the neutrino emission much weaker than that of type A or B; the transition from one 
dominating superfluidity to the other in case (CA) for vi >> 1 takes place at V2 ~ lnt>i. Moreover, 
for f 1 > 1^2 ^ 1 the factor R^^ appears to be much larger than R]^^ or R^^ (Fig- ^■ 



For not very strong superfluidities (BA) {^Jvj + v^ < 5 or [(1-T2)/0.65]''+ [(1-ti)/0.76]'' < 1), 
the reduction factor is fitted by 
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Ti)/0.8]^ < 1) Rr;j^ is fitted by the expression: 



In the case of not too strong superfluidities (CA) (v^f + i'i < 10 or [(1 - r2)/0.825]^ + [(1 
,(D) 



i?gj = 10"^ X (10" + 793.9 1;^+ 457.3 «? + 66.07 uiw? + 2.093 ?;f+ 

+ 0.3112 1;^ + 1.068^21'? + 0.01536 1;|4 + 0.006312 7;^ vf)" . (35) 

li vi = {) or V2 = 0, the above flts agree with Eqs. (|2q). The tables of -Rg^ and -R^^, as well as the 
asymptotes of these factors in the limit of strong superfluidity are given in Levenfish and Yakovlev 
(1994b). 

Before publication of the above results (Levenfish and Yakovlev, 1993, 1994b), as far as we know, 
simulation of the enhanced NS cooling with account of the superfiuidity of neutrons and protons was 
carried out in the only article by Van Riper and Lattimer (1993). The authors assumed the simplified 
reduction factor -R12 = ^1 -^2 • The actual factor, Eq. p^), is quite different, i.e., the simplified 
factor strongly overestimates the effect of superfluidity of neutrons and protons. 

As shown by Prakash et al. (1992), practically all equations of state in the inner NS cores, which 
predict appearance of hyperons, open the direct Urea reactions involving hyperons: 

Bi^B2 + l + v, B2 + l^Bi + u, (36) 

where Bi and B2 are baryons (nucleons, A-, S- or H-hyperons, A "-resonance), I is a lepton (electron 
or //-meson). The hyperonic processes ( p6| ) do not conserve strangeness; the neutrino emissivities in 
these processes are somewhat lower than in the nucleon direct Urea process (Q). The emissivities 
are given by the expressions similar to (p^), but with different numerical coefficients (Prakash et al., 
1992). 

Hyperons, like nucleons, may be in a superfluid state (Balberg and Barnea, 1998). Their super- 
fiuidity is most likely to be of singlet-state type due to their relatively small number density. The 
reduction factors for the direct Urea processes involving hyperons are similar to those for the nucleon 
direct Urea. 
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Figure 7: Lines of constant reduction factors W^' for the direct Urea in the presence of neutron and proton superflu- 
idities (AA), (BA) and (CA). The curves are labeled by the values Ig-R . In the region Tn > 1 and Tp > 1 neutrons 
and protons are normal and W^' — 1. In the region t„ < 1, Tp > 1, where only neutrons are superfluid, and in the 
region t„ > 1, Tp < 1, where only protons are superfluid, R^^' depends on the only parameter (Sect. 4.3). In the region 
T„ < 1, Tp < 1 both, neutron and protons, are superfluid at once. 
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5 Modified Urea process and neutrino bremsstrahlung due to nucleon- 
nucleon scattering 

5.1 Two branches of the modified Urea process 

In this section we will analyse the standard neutrino energy loss rates in the processes (|l|)~(^) (Sect. 
2.2). In the absence of superfluidity, these processes were studied by Bahcall and Wolf (1965a, b). 
Flowers et al. (1975), Priman and Maxwell (1979), as well as by Maxwell (1987). The latter author 
considered also the processes involving hyperons. The most detailed article seems to be that by 
Friman and Maxwell (1979), where, however, the proton branch of the modified Urea-process was 
neglected. The proton branch was analysed and the superfluid reduction factors for the standard 
process were calculated by Yakovlev and Levenfish (1995). Below in will mainly follow consideration 
of these authors. 

The main standard neutrino energy loss mechanism in a non-superfiuid matter is the modified Urea 
process. It will be labeled by upperscripts (MiV), where N = n indicates the neutron branch (|I|) of 
the process, while N = p indicates the proton branch (|2|). Both branches consist of direct and inverse 
reactions and are described by similar Feynman diagrams. If beta-equilibrium is established, the rates 
of the direct and inverse branches are equal; it is sufficient to calculate the rate of any reaction and 
double the result. The general expressions for the neutron and proton branches of the modified Urea 
process can be written as (^ = c = fee = 1): 



Q 



(MN) 



n 



\ (2vr)3 



2ee(27r)3 2e^(27r) 



e. {2^f5{Ej - E{)5{Vf - P.) tt E 1^1'' (37) 



spins 



where pj is a nucleon momentum (j = 1, 2, 3, 4), pe and £e are, respectively, the electron momentum 
and energy; while py and £i, are, respectively, the neutrino momentum and energy. The delta function 
5{Ef—Ei) describes energy conservation, and 5(P/— Pj) describes momentum conservation; subscripts 
i and / refer to the initial and final particle states, respectively. Furthermore, C means product of 
the Fermi-Dirac functions or corresponding blocking factors for nucleons and electron; \M\^ is the 
squared reaction amplitude, the sum is over all particle spins. The factor 2 in the denominator before 
the summation sign is introduced to avoid double counting of the same collisions of identical particles. 
Introducing dimensionless variables, the emissivity (|37| ) can be written (Yakovlev and Levenfish, 
1995) in the form similar to the direct Urea emissivity (see Eq. (pO|) and explanations afterwards): 
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The quantity A contains integrals over orientations of particle momenta (j=5 corresponds to electron); 
all lengths of the momenta pj of nucleons and electron in the delta function should be set equal to the 
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Fermi momenta. Typical neutrino momentum pi, ~ T is much smaller than the momenta p-p of other 
particles. Thus, it can be neglected in momentum conservation; the integration over orientations of the 
neutrino momentum in A immediately yields Air. The quantity /, given by Eq. (^0|), includes integrals 
over dimensionless energies of neutrino x,^ = Pu/T = SyjT and other particles Xj = v^Ap — pp-)/T 
(see Eq. (|8|)); fj = [exp(xj) + 1]~^. Finally, the quantity S is the product of density of states of the 
particle species 1 < j < 5 (with effective masses m*,) at the Fermi surfaces. 

Owing to the similarity of the Feynman diagrams for the neutron and proton branches of the 
modified Urea process, ([l|) and (P), the reaction amplitudes M and phase integrals / are practically 
the same. The reactions differ by the products of the density of states S and by the angular integrals 
A. 

In a non-super fluid matter, one has I = Iq = 11513 vr^ /120960 (see, e.g., Shapiro and Teukolsky, 
1983). The amplitude M for the neutron branch of the modified Urea process was calculated by 
Friman and Maxwell (1979) using the Weinberg-Salam-Glashaw theory of electroweak interaction. 
Long-range (small momentum transfer) part of the nucleon-nucleon interaction was described in the 
one-pion-exchange approach, while the short-range (large momentum transfer) part was described in 
the frame of the Landau Fermi-liquid theory (e.g., Baym and Pethick, 1991). The squared amplitude 
|Mp summed over spin states is given by Eq. (39) in Friman and Maxwell (1979). 

The angular integrals (|3^) for the processes ([l|) and (^) differ due to the difference of neutron and 
proton Fermi momenta (Sect. 2.1). The proton reaction branch involves three protons with rather small 
Fermi momentum and the only neutron with larger momentum, while the neutron branch involves 
three neutrons with large momentum and the only proton with small momentum. Direct calculation of 
the angular integrals for these processes in the absence of superfluidity gives (Shapiro and Teukolsky, 
1983; Yakovlev and Levenfish, 1995): 






^no = ^^T^, ^pO = ^V^(l-^le, (42) 



where = 1 if the proton branch is allowed by momentum conservation, and = otherwise. Let 
us remind that in the outer NS core the Fermi momenta pp^ and p-p^ are much smaller than pp^ (Sect. 
2.1). Thus O = 1 for pp„ < 3pp +Pp(.- Notice that the expression for Ano is obtained from Eq. ( [39| ) in 
case pp„ > Ppg + Ppp- In the opposite case oi pp^ < Pp^+Ppp, Ano is given by Eq. (13) in Yakovlev and 
Levenfish (1995), but in this case the direct Urea process dominates and the modified Urea processes 
are insignificant. 

Using the above expressions for Iq and Ano, Friman and Maxwell (1979) calculated the neutrino 
emissivity in the neutron branch of the modified Urea process (in the standard physical units) : 

(Mn) _ 11513 Glglmfm; ( f- ^ PFei^BTf 
^^ ~ 30240 2^ \mj h'^'c^ ""^" 



.55 X 10^^ ( :::^ ) ( i:i! )( iifi ) '\« a./?, ^i^, (43) 

\mnJ \mpj \noJ " cm^-s 

where gA = 1-26 is the axial-vector constant of weak hadron current, /'^ ~ 1 is the vrA^-interaction 
constant in the p-state in the one-pion-exchange approximation, and ttt-tt is the pion mass (tt ). The 
factor a„ describes momentum dependence of the squared reaction amplitude in the Born approxi- 
mation, and Pn contains different corrections. According to Eq. (62) in Friman and Maxwell (1979) 
a„ ~ 1.76 — 0.63 (rao/n„)^'^, where n„ is the number density of neutrons. In their final Eq. (65c) for 
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Qq " Priman and Maxwell (1979) used the value a„ = 1.13, calculated for p = po, and set /?„ = 0.68 
(to account for the correlation effects). 

The expression for the neutrino emissivity in the proton reaction was derived by Yakovlev and 
Levenfish (1995): 
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(44) 



Taking into account all uncertainties concerned with the calculation of the reaction amplitude by 
Friman and Maxwell (1979) we set Op = a„ and Pp = Pn- It is easy to see that the emissivities in the 
neutron ( |43| ) and proton ( |4^ ) branches of the process are similar. The main difference of the proton 
branch is in its threshold character: the reaction is allowed for pp^ < 3pp + pp^. In npe-matter, this 
inequality is equivalent to pp„ < 4ppg, i.e., to He > Un/^^- The latter condition is realized almost 
everywhere in the NS core. It can be violated only for ultra-soft equations of state al p < pQ. For 
these equations of state, the proton branch can be forbidden in the outermost part of the NS core. 
Similar but much more stringent threshold conditions can be formulated for the direct Urea process 
(D (Sects. 2.2 and 4). 

Comparing Qq ^ and Qq " , we find: Qq ^ /Qq " = (m*/m,*)^[l —Ppf./{4:Pp)]. For instance, at 

m* = m* and Ppg= p-pp we have Qg ^ = 0.75 Qq " • Therefore the proton branch of the process is 
nearly as efficient as the neutron branch. 

The potential efficiency of the proton branch was outlined by Itoh and Tsuneto (1972) who, 



,(Mp) 



(Mp) 



however, did not calculate Qq . Later the neutrino emissivity Qq was calculated by Maxwell 
(1987) who found it negligibly small as compared to Qg " . This conclusion is erroneous due to 
several inaccuracies made by Maxwell (1987) (and analysed by Yakovlev and Levenfish, 1995). In 
particular. Maxwell incorrectly neglected the electron momentum in momentum conservation. 



5.2 Modified Urea proeesses in superfiuid matter 

Nucleon superfluidity reduces the modified Urea processes (see Sect. 4.2). Let us analyse this reduction. 
We will adopt the traditional assumption (Sect. 3.1) that the proton superfiuidity is of type A, while 
the neutron superfiuidity is of type B. Our results will also be valid for the neutron pairing of type 
-^ (^-t /5 < Po)) as will be mentioned later. Reduction of the modified Urea by the neutron pairing of 
type C has not been considered so far. 



Superfiuidity affects the nucleon dispersion relations under the integral (37) in accordance with 
Eq. (^ . The neutrino emissivity of the neutron and proton branches of the process can be presented 

in the form 

^(Mn) ^ g(Mn)^{M„)^ ^(Mp) ^ gl-P;^U-P; ^ (45) 

and ( ^4[ ) in non-superfiuid matter, while W "' and 
R} ^' are the factors, which describe superfiuid reduction of the reactions {w < 1). Generally, 
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where Qq " and Qq ^ are the emissivities (li^) 
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(46) 
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f{z) = l/(e^ + 1); Zj for nucleons {j < 4) is defined by Eq. (|^; one should set 25 = X5 = Xe for an 
electron. 

Equation ( |4^ ) enables one to calculate the reduction factors w "' and R^ ^' as a function of T, 
Ten and Tcp. Below we present the results (Yakovlev and Levenfish, 1995) for the proton superfluidity 
of type A and normal neutrons as well as for the neutron superfluidity of type B and normal protons. 
The behaviour of W "' and R^ ^' in the presence of the joint superfluidity of nucleons is described 
in Sect. 5.4. 



Singlet-state proton pairing. Since the singlet-state gap is isotropic the angular and energy 
integrations in Eq. (^6|) are decomposed, and the angular integration remains the same as in non- 
superfluid matter. 

Just as in the case of the direct Urea, the reduction factors of the neutron and proton branches of 
the modified Urea process, -R!,^" and -Rp^^^ , can be expressed in terms of the dimensionless energy 
gap Vp = VA (see Eq. (|l^); here and hereafter the subscripts in R^ ' indicate the superfluid particle 
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species and superfluidity type). Clearly, i?^^ = 1 and Rpj^ =1 for T > T^p (u^ = 0). In the case 
of strong superfluidity (T <^ Tcp, v^ -^ 00), the asymptotes of the both factors are (Yakovlev and 
Levenfish, 1995) 
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(47) 
(48) 
(49) 



where Tp = T/Tcp. 
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Yakovlev and Levenfish (1995) calculated i?!,^ and R^ for intermediate v^ and fitted them by 
the analytic expressions, which reproduced also the asymptotes in the limit of f ^ 00 and obeyed the 
condition ii(^^)(0) = l: 
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Equations ( |50| ) and (^), together with (jT^), fully determine the dependence of -Ri^" and Rp^^ on 



The above results (Yakovlev and Levenfish, 1995) are valid also for the singlet-state superfluidity 
of neutrons. Evidently, in that case one should set f„ = va and 



R^ 



(Mp), 



?(Mn) 



(Mn), 



JMp), 



tT'ivA) = R';r{vA), KTivA) = R';rivA) 



(52) 



Then Eq. (50) describes reduction of the proton branch of the process and Eq. (51) describes reduction 
of the neutron branch. 
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Wolf (1966) (Fig. 2 of his article) as well as Itoh and Tsuneto (1972) were the first who considered 
the reduction factor -R„^" of the neutron branch of the modified Urea process by the singlet-state 
neutron superfluidity; this factor is analogous to the factor R^j^^ ■ Notice that Itoh and Tsuneto 
(1972) analysed only the asymptote (|48|). In both articles the same asymptote (48) was obtained but 
with different numerical factors ^. Wolf (1966) got ^ = 0.123, while Itoh and Tsuneto (1972) obtained 
^ = 7r/15 ss 0.209. Recently -R^^" has been calculated independently by Pizzochero (1998) under 
artificial assumption that the superfluid gap is temperature independent. His results are described by 



the factor ^„^" (^_a) presented above with Vy^ = 1.764/t. 



Triplet-state neutron pairing. In this case the neutron gap is anisotropic. The proton branch 
of the modified Urea process is analysed easily since the only one superfluid particle is involved. The 
expression for the reduction factor -R„^ reduces to a one-dimensional integral over angle ■!?„ (between 
neutron momentum and quantization axis) of the factor RS^'^ (v) fitted by the expression (pOJ). An 



'-pA 
argument va in the latter expression should be formally replaced by ys in accordance with Eq. (ITOl) 

(Sect. 3.2). It is evident that R^b \^) = 1 for u = % = 0. In the limit % ^ oo one can use the 

asymptote ( p7[ ) in the integrand. Then for T <C Ten (% ^ 1) according to Yakovlev and Levenfish 

(1995) 

„{Mp) 72 7 3.99 X 10-6 / 1.188\ 

""-' = 11513 vr^Vs"^'"^^"'^^ = rl ^"^ l"^ ' ^''^ 

where r„ = T/Tcn- Yakovlev and Levenfish (1995) calculated also R^^ for intermediate values of % 
and fitted the results by the analytic expression: 

,7 I f,5 



^(M,) ^ a_^ ^^p 1^2 3gg _ ^(2.398)2+^2^ 



(54) 



0.1612 + J(0.8388)2 + (0.1117t;5)2, b = 0.1612 + J(0.8388)2 + {Q.UI^VbY 



Exact calculation of the reduction factor -R„^" of the neutron branch of the modified Urea process 
by the triplet-state neutron superfluidity for intermediate values of % is complicated; an approximate 
expression will be given in Sect. 5.4. Here we present only the asymptote of R^^^ in the limit r„ ^ 1 



(Yakovlev and Levenfish, 1995) 

^"^"^ = 11^ X ^^v-s-M-2v,) = —— exp [-^-^). (55) 



(Mn) 120960 2 . « , „ . 1.56 x 10"^ / 2.376 



In this case, as for the proton reaction involving superfluid protons (48), the effect of superfluidity 
appears to be very strong: the exponent argument in -R^^" contains the doubled gap (three neutrons 
which participate in the reaction belong to superfluid component of matter). 

The dependence of the reduction factors Rp^^ , RpA^ ^nB^ ^^^ ^nB^ "-"^ dimensionless tem- 
perature T and dimensionless gap parameter v is plotted in Figs. |8l and ^. For comparison, we present 
also the reduction factors of the direct Urea process, -R^ and Rj^ . Let us mention that in the case 
of the strong neutron superfluidity (T <^ Ten) and normal protons the proton branch of the modified 
Urea process becomes much more efficient than the neutron branch. However in this case the main 
neutrino emission comes from the neutrino bremsstrahlung due to pp-scattering, which is not affected 
by the neutron superfluidity. 
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Figure 8: Reduction factors of various neutrino emission processes by superfluidity of neutrons or protons versus T/T^. 
Curves 1 show reduction of the p-branch of the modified Urea process (solid line) and the direct Urea process (dashed 
line) by neutron superfluidity of type B. Curves 2 correspond to reduction of the n- branch of the modified Urea (solid 
line) and the direct Urea (dashed line) by proton superfluidity of type A. Dot-and-dash lines 5, 6 and the solid line 3 
refer to the np, pp-scattering and p-branch of the modifled Urea processes, respectively, for the same superfluidity. Solid 
line 4 is the asymptote of the reduction factor for the n-branch of the modified Urea process due to n superfiuidity of 
type B. 

5.3 Neutrino bremsstrahlung due to nucleon-nucleon scattering in superfluid mat- 
ter 

Let us consider the effect of superfluidity on generation of neutrino pairs in nn, np and pp-scattering @ ; 
these processes are often referred to as neutrino nucleon-nucleon bremsstrahlung. In a non-superfluid 
matter, the neutrino bremsstrahlung is about two orders of magnitude less efficient than the modified 
Urea process. Corresponding neutrino emissivities (in the standard physical units) are (Friman and 
Maxwell, 1979; Yakovlev and Levenfish, 1995) 
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Figure 9: Same as in Fig. pi but versus dimensionless gap parameter. 
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(58) 



where m-j^ is the n^ mass, and Mu is the number of neutrino flavors. aj\f]\[ and Pnn have the same 



meaning as in Eqs. (43) and (|44|); they are slowly varying functions of nucleon Fermi momenta, i.e., 
density. Friman and Maxwell (1979) assumed the factors unn to be equal to their values at p = p^: 
1.06, Upp = 0.11, and the correction factors (line d = 0.7 fm in their Table 1) 



ar. 



0.59, a 



np 



Pnn to be Pnn = 0.56, Pnp = 0.66. They did not calculate /3„„. The presented values of /?„„ and p,, 



-'np 



indicate that it is reasonably to set Ppp ~ 0.7. Hereafter we will use Mu = 3 in Eqs. (|5q)-(p8|) while 
Friman and Maxwell (1979) took into account two neutrino flavors. 

In analogy with the Urea processes let us introduce the factors R^^'^), which describe superfluid 
reduction of the neutrino bremsstrahlung due to A^A^ scattering: 



q{NN) ^ q{NN)j^(nN) 



(59) 



Consider a singlet-state pairing A of neutrons or protons. Then the modification of the dispersion 
relation by superfluidity affects only the phase integral. For a non-superfluid matter, the latter integral 
is given by Eq. (46) in Friman and Maxwell (1979) (cf. with Eq. (pQ)): 
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where in this case Xj^ = Ey/T determines the total energy £y of a neutrino pair. Accordingly the 
reduction factors take the form: 



^ ~ 1647r8 io ' " 



n 



+00 



dxj f{zj) 



ME 
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where the dimensionless quantity Zj is defined in Eq. (^ . It is clear that R^^"> = 1 for r > 1 . 

In the case of the singlet-state proton pairing the reduction factors for the np and pp processes, 
Rp^ and Rp^ , are reduced to two dimensional integrals which can be calculated numerically 
(Yakovlev and Levenfish, 1995). In the limit of strong superfluidity (rp ^ 1, u^ ^ 00) the asymptotes 
of these factors are: 
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where ^1 ss 849. The asymptotes (|6^) and (|63|), as well as numeric values of R^j^ and R^X calculated 
for intermediate values of f^ can be fitted by the expressions: 
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a = 0.9982 + J(0.0018)2 + (0.3815^^)2, b = 0.3949 + a/(0.6051)2 + (0.2666t;^)2 ; 
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For the singlet-state neutron pairing, we evidently have 
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The reduction factors considered above are shown in Fig. ^ versus T/T^ and in Fig. |9| versus 
dimensionless gap parameter v. 

Summarizing the results of Sects. 5.2 and 5.3 we note that the exponent argument in the asymp- 
tote of a reduction factor contains a single gap if one or two reacting particles belong to superfluid 
component of matter. Under the number of reacting particles we mean total number of particles (in 
the initial and final states of a bremsstrahlung process, in the initial and final state of direct or inverse 
reaction of an Urea process) belonging to superfluid component. The gap in the exponent argument 
is doubled if three or four superfluid particles are involved, etc. As seen from Fig. ^, the factor Rl^ 
(two superfluid particles) falls down with increasing the superfluidity strength (with increasing v) 
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much more rapidly than the factors -R!,^" or R\^^ (one superfluid particle). Accordingly, R!^^ (four 

superfluid particles) falls down faster than Rpj^ (three superfluid particles). Let us stress that we 
discuss the neutrino reactions considered in Sects. 4 and 5. The formulated rule is invalid for the 
neutrino emission due to Cooper pairing of nucleons (Sect. 6.1). 

5.4 Neutrino reactions in the presence of neutron and proton superfluidity 



If neutrons and protons are superfluid at once calculations of multi-dimensional integrals ( |4q ) and (plj), 
which determine reduction of the standard neutrino reactions, becomes very complicated. However, if 
a very high accuracy is not required, one can avoid calculation by noticing that the reduction factors 
for the processes involving one superfluid particle {R\ , -Ri^" , Rg and R\^b ) ^'^^ close to one 
another as functions of the dimensionless parameter v (Fig. P). This enabled Levenfish and Yakovlev 
(1996) to formulate approximate similarity criteria for different reduction factors. Using these criteria 
we have constructed the approximate reduction factors for the proton and neutron branches of the 
modified Urea process in the presence of neutron and proton superfluidity: 

Wba {Vn,Vp) W (D), , RuB i^n), (67) 

Rb y"n) 



R^li'^Vn^Vp) 



RbA {Vn,Vp) PS (D) RpA yVp). (68) 

Ra W 

We expect that these factors, as functions of corresponding parameters v (corrected due to the number 
of superfluid particles) do not differ strongly from the reduction factor i?^^(t'„, Vp) for the direct Urea 
process. If protons are normal {vp = 0), then the expression for R^J' {vn,Vp) becomes exact; if neutron 
are normal {vn = 0), the factor Rg^ {vmVp) is exact. In addition, the approximate factors satisfy 



the relationship analogous to (l33|). 

One can also expect similarity of the reduction factors for other neutrino reactions. For instance, 
the reduction of the neutron branch of the modified Urea process, R^b^ , by a moderate neutron 
superfiuidity {v < 10) should not deviate strongly from the reduction of the proton branch of the 

modified Urea, Rpj^^ , by the proton superfluidity: 

4r-<^^K). (69) 

The approximate reduction factor R^ for the nn scattering in the presence of n superfluidity and 
the approximate reduction factor Rj^X ^^ the np scattering by the neutron and proton superfluidities 
can be written as 

Ri7^ - R^^fM, (70) 

p(D), ^ 

RbA - ID), . RpA i^pi- {'^) 

Ra W 



In the absence of the neutron superfluidity, Eq. (|7lD becomes exact. 

The modified Urea and neutrino bremsstrahlung processes can involve hyperons. After minor 
modifications, the above reduction factors can be valid for these reactions as well. 
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6 Neutrino emission due to Cooper pairing of nucleons 

6.1 Neutrino emissivity produced by nucleon pairing 

In contrast to the neutrino emission processes considered in Sects. 4 and 5, this process is ahowed only 
in the presence of superfluidity (Sect. 2.2): the superfluidity distorts the nucleon dispersion relation 
near the Fermi surface and opens the reaction (|5|). Actually the process consists in emission of a 
neutrino pair by a nucleon whose dispersion relation contains an energy gap; however, in theoretical 
studies, it is convenient to use the formalism of quasi-particles and treat it (Flowers et al., 1976) as 
annihilation of two quasi-nucleons N into a neutrino pair: 

N + N ^u + i?. (72) 

The reaction goes via weak neutral currents and produces neutrinos of all flavors. Following Yakovlev 
et al. (1999) we will outline derivation of the neutrino emissivity due to singlet-state or triplet-state 
pairing of non-relativistic nucleons. The reaction is described by the Hamiltonian (h = c = k^ = ^) 

H = -^{cvJoIo-caJI), (73) 

where Gp is the Fermi constant, and cy and ca are, respectively, the vector and axial-vector constants 
of neutral hadron currents. For neutron currents, we have (see, e.g., Okun', 1990) cy = 1, ca = gA = 
1.26, while for proton currents cy = 4 sin^ 0-w — 1 « —0.08, ca = —gA, where ©w is the Weinberg 
angle, sin ©w = 0.23. Strong difference of cy for neutrons and protons comes from different quark 
structure of these particles. Furthermore, 

J'^ = ( jO, j) = (^+^, 4>+a^) , I'' = V^,7'^(l + 7^)^;, (74) 



are 4-vectors of neutral currents of quasi-nucleons and neutrino (/i=0,l,2,3), respectively; ipu is the 
neutrino wave function, upper bar denotes Dirac conjugate; 7^ and 7^ are Dirac gamma-matrices, a is 
a vector Pauli matrix; ^ is a second-quantized quasi-nucleon wave function. The function ^ is derived 
using the Bogoliubov transformation. Its description for the singlet-state and triplet-state pairing is 
given, for instance, in Lifshitz and Pitaevskii (1980) and Tamagaki (1970). In the both cases 

^ = E :^- [e-"*+'P'-C/<.,(p)ap, + e-*-'P'-V.,(-p)a+J , (75) 

pcrr; 



where p and e = \/vp{p — Pp)'^ + 6p are, respectively, the quasiparticle momentum and energy (with 
respect to the Fermi level). A basic spinor Xa describes a nucleon state with fixed spin projection 
(o" = ibl) onto the quantization axis (axis z); r] enumerates quasi-nucleon spin states; Sp is the energy 
gap in the quasi-particle spectrum, Vp is the Fermi velocity, a^^ and cip^ are, respectively, creation 
and annihilation operators for a quasiparticle in a (prj) state; U{p) and V^(p) are the operators of 
the Bogoliubov transformation. For |p — Pp| *C Pp, their matrix elements obey the relationships 
U(Tr]{p) = UpSar/ and J2ari l^ar?(p)|^ = 2 v^, where 
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(76) 



For a singlet-state pairing, the gap 5p is isotropic, and the quantities Up and Vp depend only onp= |p| 
For a triplet-state pairing, the quantities 5p, Up and Vp depend also on orientation of p. 
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Let Qi, = (cjjy, qjy) and ql = (col, q'j^) be 4-momenta of neutrino and antineutrino, respectively, while 
p = (e, p) and p' = (e', p') be 4-monienta of annihilating quasi-nucleons. Using the Golden Rule of 
quantum mechanics, we can present the neutrino emissivity due to Cooper pairing (CP) as: 



Q^'''' = (^) l^J^M^mm 



\2V2J 2'^" J (2^)3 (2^)3 

f dq^ dq[, r , ,, ,2 , 2 . , , 
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(77) 



where Afiy=3 is the number of neutrino flavors, and the factor 1/2 before Mu excludes double counting 
of the same quasi-nucleon collisions. The integral is taken over the range {q,^ + q[,)'^ > 0, where the 
process is open kinematically; /(e) = l/[exp(e/T) + 1], i,k = 1,2,3; 

/oo = E K^I^+'I'I^)!', hk = E (^l^+^.'I'l^) {Bl^^akMA)*. (78) 

rjr]' r)ri' 

Here 1^4) stands for an initial state of the quasi-particle system in which one-particle states (p, rj) and 
(p', r]') are occupied, and \B) stands for a final state of the system in which the indicated one-particle 
states are empty. 

The integral ( |77| ) is simplified by the standard technique, as described in Yakovlev et al. (1999). 
The transformations take into account that nucleons are non-relativistic and strongly degenerate as 
well as the fact that the process is open kinematically in a small domain of phase space where the 
quasi-nucleon momenta p and p' are almost parallel. The latter circumstance allows one to take 
smooth functions Ioo(p,p') and /jjt(p,p') out of the integral over dp' putting p' = —p. After a 
number of transformations the final expression for the neutrino emissivity can be written (in the 
standard physical units) as (Yakovlev et al., 1998, 1999) 

= 1.170 X 10^1 f^) (^^\ TlM, aF ^^, (79) 

yrriN J XniNcJ curs 

where Tg = T/(10^ K), a is a numeric factor (see below), and the function F, in our standard notations 
(|8|), is given by the integral 

F = ^ /dOy2 r^!^. (80) 

47r i ^ io (e^ + 1)2 ^ ^ 

The singlet-state gap 6p is isotropic; thus, integration over dil is trivial and gives 47r. In the triplet- 
state case, the function F contains averaging over positions of a quasi-nucleon on the Fermi-surface. 
While using Eq. ([79|), one can take into account that p-p/i'mj^c) « 0.353 {n^ /rioY'^ , where tin is the 
number density of nucleons N , uq = 0.16 fm~3_ 

The emissivity Q depends on superfluidity type through the factor a and the function F. For 
the singlet-state neutron pairing, a is determined by the only vector constant cy: anA = Cy = ^- If 
we used similar expression for the singlet-state pairing of protons, we would obtain a very small factor 
flpA = 0.0064, which indicates the weakness of the process. Under these conditions, one should take 
into account the relativistic correction to a, produced by the axial-vector proton current. Calculating 
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and adding this correction for the singlet-state pairing of protons, Kaminker et al. (1999b) obtained 



apA = cy + c\ 



^Fp 



m^ 



m„ 



+ 
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0.0064 + 1.59 ' ^^ 



m„ 



nir. 



+ 
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(81) 



where v-p Jc = {pp /mpc){mp/m*). The relativistic correction appears to be about 10 - 50 times 
larger than the main non-relativistic term. This enhances noticeably the neutrino emission due to the 
singlet-state proton pairing although it remains much weaker than the emission due to the neutron 
pairing. 

In the case of the triplet-state pairing a is determined by the both, vector and axial-vector, con- 
stants of neutral hadron currents: a = a^B = clnc = Cy + 2c\ (Yakovlev et al., 1999). For the neutron 
pairing, we obtain UnB = o-nC = 4.17. Notice that in the case of the triplet-state pairing of protons 
which is thought to be hardly possible in the NS cores we would obtain Ups = ipC = 3.18. Under 
such exotic conditions, the neutrino emission due to proton pairing would be almost as efficient as the 
emission due to neutron pairing. 

The result for the singlet-state pairing of neutrons, presented above, coincides with that obtained 
in the pioneering article by Flowers et al. (1976) (for two neutrino flavors, Mu = 2). Similar expressions 
obtained by Voskresensky and Senatorov (1986, 1987) for Mu = 1 contain an extra factor (1 + 3g\). 
In addition, the expression for Q^ ' derived by the latter authors contains a misprint: vr^ in the 
denominator instead of vr^ (although numeric estimate of Q^ ' is obtained with the correct factor 
vr^). The cases of the singlet-state proton pairing and the triplet-state neutron pairing have been 
analysed by Yakovlev et al. (1998, 1999) and Kaminker et al. (1999b) for the first time. 

Let us mention that the values of cy, ca and a can be renormalized in dense NS matter under 
in-medium effects. The renormalization is a difficult task which will be neglected below. 

The function F, given by Eq. (|80|), depends on the only variable v, the gap parameter (see Eq. 
@). Using Eq. ( pO| ) one can easily obtain the asymptote of this function and calculate its dependence 
on r = T/Tc for superfluidity of types A, B and C in analogy with calculations presented in Sects. 
3-5. This has been done by Yakovlev et al. (1998, 1999). 

In a small vicinity of T ~ Tc, in which w <C 1 and r — > 1, we have: 



Fa{v) = 0.602 1;2 =5.65(1 -r), 
Fb{v) = 1.204 f 2 = 4.71 (1 - r), 
F(j{v) = 0.4013 1;2 = 4.71 (1 - r) 



(82) 



For the low temperatures T <^Tc the parameter v ^ 1, and the asymptotes of F{v) are 



Fa{v) = ^ v'^'^l'' exp(-2z;) = -^^ exp 



^ 13/2 



Fn(v) 



Fdv) 



vr g 1.27 

V exp(— 2f ) = — 7- exp 



4^3 
50.03 



2.376 



(83) 



12. Ir^ 



Let us stress that the neutrino emission due to nucleon pairing differs from other neutrino reac- 
tions: first, is has a temperature threshold (is allowed only for T < Tc); second, its emissivity is a 
nonmonotonic function of temperature. The emissivity grows rapidly with decreasing T just after 
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Figure 10: Temperature dependence of the neutrino emissivity due to Cooper pairing of neutrons for p = 2 x 10^* g 
cm~^ and Ten = lO" K for superfluidity types A (solid line), B (dashes) and C (dots). 

superfluidity but then reaches maximum and decreases. According to Eq. (|83|), a strong superfluidity 
reduces considerably the emissivity just as it reduces the heat capacity or direct Urea process: the 
reduction is exponential if the gap is nodeless (cases A and B) and it is power-law otherwise (case 
C); see Sects. 3 and 4. 

The asymptotes ( |8^ ) and (p^), as well as numeric values of F(y) for intermediate u, are fitted by 
the simple expressions (Yakovlev et al., 1998, 1999) 

1/2 



F^{v) = (0.602 u^ + 0.5942 w"^ + 0.288 v*') 0.5547 + J (0.4453)2 + 0.0113^2 



Fn(v 



X exp 1 --^4^2 + (2.245)2 + 2.245 

1.204 1;2 + 3.733 v^ + 0.3191 v^ 
1 + 0.3511 7;2 



0.7591 + J (0.2409)2 + 0.3145^2 



(84) 



Fdv) 



exp ( -y^4t;2 + (0.4616)2 + 0.4616 j , 

0.4013 ^2 - 0.043 u^ + 0.002172 v^ 
' 1 - 0.2018 ^2 + 0.02601 v^ - 0.001477 v^ + 0.0000434 ', 



Equations (79) and (Q) enable one to calculate easily the neutrino emissivity Q^ ' due to Cooper 
pairing of nucleons for superfluidity of types A, B and C. Similar expressions describe neutrino 
emissivity due to pairing of hyperons. The required values of a are listed in Yakovlev et al. (1999). 

Figure |l^ (from Yakovlev et al., 1999) shows temperature dependence of the emissivity Q^ ' due 
to neutron pairing in the NS core for p = 2 x 10^^ g cm~^. The adopted equation of state of matter 
is described in Sect. 7.2. The effective nucleon masses are set equal to m^ = 0.7m]\[, and the critical 



41 



o 



bJO 



bJO 





_ 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 


' ' ' . 


21 


- T,„=4x10'k 


^ 


^ ^ 




" T,p=6xlO^K 




^ 


20 




^ 


^ 


^^ 


19 


— - ■ ' ' / / ^^/ 




- 




- /^ /n 




- 


18 


- /p / 




- 




" i/ / , , 1 





- 



8.5 



9.0 

IgT [K] 



9.5 



Figure 11: Temperature dependence of the total and partial (n, p, e) specific heat capacities at p = 2po for neutron 
superfluidity of type B with Ten = 4 x 10^ K and proton superfluidity of type A with Tcp = 6 x 10* K. Dotted lines 
show corresponding heat capacities in non-superfluid matter. 

temperature is Ten = 10^ K. The density of study is typical for transition between the single-state 
pairing and the triple-state one (Sect. 3.1). Thus, different models of nucleon-nucleon interaction may 
lead to different neutron superfluidity type. We present the curves for the"three superfluidity types 
considered above. 

When the temperature falls down below T^n the neutrino emissivity produced by Cooper pair- 
ing strongly increases. The main neutrino energy release takes place in the temperature interval 
0.2 Ten < r < 0.96 Ten, with the maximum at r~0.4Tcn. The emissivity may be sufficiently high, 
compared with or even larger than the emissivity of the modified Urea process in non-superfluid matter 
(Sect. 4). Under certain conditions, neutrino emission due to pairing of neutrons may be significant in 
the inner NS crust (Yakovlev et al., 1998; Kaminker et al., 1999). The reaction may be noticeable even 
in the presence of the direct Urea process in the inner NS core if the direct Urea is partly suppressed 
by the proton superfluidity (see below). 



6.2 Summary of Sections 3 — 6 

Let us summarize the results of Sects. 3-6. For illustration, we use the same equation of state of 
matter (neutrons, protons and electrons) in the NS core as in our cooling simulations (Sect. 7.2), and 
set m^ = 0.7 rriAr. Let us adopt that the neutron pairing is of type B, while the proton pairing is of 
type A. 



Figure 11 illustrates the effect of superfluidity on the heat capacity for p = 2pQ, Ten = 4 x 10^ K 



and Tcp = 6 x 10^ K. In the absence of superfluidity, the main contribution into the heat capacity 
comes from neutrons. The heat capacities of protons and electrons are lower than that of neutrons 
by a factor of 2.5 and 9, respectively. After superfluidity appears with decreasing temperature, the 



42 




S''' 



9.0 9.5 

IgT [K] 



8.5 9.0 

IgT [Kl 



Figure 12: Temperature dependence of the neutrino emissivity in different reactions for the neutron superfluidity of 
type B with Tc„ = 8 x 10^ K and the proton superfluidity of type A with Tcp = 4 x 10^ K at p — 2 po (standard 
neutrino reactions, Fig. a) and p = 5 po (direct Urea is aUowed, Fig. b). Dot-and-dash line shows the emissivity due to 
Cooper pairing of neutrons plus protons; solid line presents the total emissivity. In Fig. a the dashed line gives the total 
emissivity in two branches of the modified Urea process; the dotted line exhibits the total bremsstrahlung emissivity due 
to nn, np and pp scattering. Figure b: the dashed line corresponds to the direct Urea process. 

relative contributions of different particles change. The jumps of the total and neutron heat capacities 
with the fall of temperature below T = Ten are associated with latent heat release produced by pairing 
of neutrons. However for T < 10^ K the neutron superfluidity becomes strong and reduces greatly the 
partial heat capacity. The main contribution into the heat capacity comes now from protons. When 
the temperature falls below Tcp, the proton heat capacity up jumps due to appearance of the proton 
superfluidity. For T < 2 x 10*^ K the latter superfluidity, in its turn, becomes strong and reduces 
exponentially the proton heat capacity. As a result, at lower temperatures the total heat capacity 
is determined by electrons and does not depend on nucleon superfluidity. Examining Fig. |ll|, it is 
easy to predict relative contributions of different particles into the heat capacity for any relationships 
between T, Ten and Tcp. 



The effect of superfluidity on neutrino reactions is more complicated. For instance, Figs. 12 show 



the neutrino emissivities in different reactions for Ten = 8 x 10 K and Tcp = 4 x 10 K. Figure 



12a corresponds to p = 2pQ. The direct Urea process is forbidden at this density (being allowed at 
Per = 4.64 Po = 1.30 X 10^^ g cm^^, for a given equation of state). In the absence of the neutron 
superfluidity (T > Ten) the dominant mechanism is the modified Urea process. If, however, the 
temperature decreases from T = Ten to T ~ 10^'^ K the total neutrino emissivity increases by about 
two orders of magnitude due to Cooper pairing of neutrons. Therefore, sometimes the appearance of 
superfluidity accelerates NS cooling instead of slowing it. 

Figure \l^ corresponds to a denser matter, p = Spoi where the powerful direct Urea process is 
switched on. In this case, the neutrino emissivity is actually determined by two processes, the direct 
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Figure 13: Regions of Ten (superfluidity of type B) and Tcp (type A), in which different neutrino reactions dominate 
at T = 10®, 3 X 10* and 10* K in matter of density p — 2 po (standard cooling) and p — 5 po (enhanced cooling). 



Urea and Cooper pairing of nucleons. The direct Urea dominates at T > 3 x 10® K. With further 
decrease of T, the direct Urea and the reactions considered in Sects. 4 and 5 are reduced so strongly, 
that the Cooper pairing becomes dominant. 

It is well known that the neutrino emission from a non-superfluid NS core is mainly determined 
by a single, most powerful neutrino emission mechanism: the direct Urea for the enhanced cooling, 



or the modified Urea for the standard cooling. However, as seen from Fig. 12, this "simplicity" is 
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violated in the superfluid NS cores. Different neutrino mechanisms can dominate at different coohng 
stages depending on T, Ten, T^p, and p. 

Figures ^ show which neutrino mechanisms dominate for different Ten and T^p- In addition to 
the neutrino processes considered above we have included one more process, neutrino bremsstrahlung 
due to electron-electron collisions (Kaminker and Haensel, 1999). This mechanism, as a rule, is 



sufficiently weak and neglected in cooling simulations. Three left Figs. |1J illustrate the case of the 
standard cooling at p = 2pQ for three internal stellar temperatures: 10^, 3 x 10*^ and 10^ K; three right 
figures correspond to the enhanced cooling at p = 5 po for the same T. The chosen values of T cover 
the region most interesting for practice. Our calculations show that topology of the figures varies only 
slightly with p as long as p does not cross the threshold value p = pcr- Therefore, the presented figures 
reflect adequately the efficiency of all neutrino processes in the core of the cooling NS. One can see 
that, in the presence of superfluidity, many different mechanisms can dominate in certain parameter 
ranges. 

Notice that if the direct Urea is open, the modified Urea is always insignificant and may be ne- 
glected. In the presence of the neutron superfluidity alone, the neutrino bremsstrahlung due to pp 
collisions becomes the main mechanism at T <C Ten being independent of the neutron superfluidity. 
In the presence of the proton superfluidity alone, the neutrino bremsstrahlung due to nn collisions 
dominates at T ^ Tcp. The Cooper pairing of neutrons exceeds the standard neutrino energy losses 
for T < 10^ K and for not too strong neutron superfluidity (0.12 < T/T^n ^ 0.96). This parameter 
range is very interesting for applications. The neutrino emission due to Cooper pairing of neutrons is 
also significant at early cooling stages, when T > 10^ K, but in a narrower temperature range near 
T~ 0.4 Ten, or in the presence of the proton superfluidity. Although the neutrino production due to 
proton pairing is much weaker, it can also dominate. The neutrino emission due to pairing of neutrons 
and protons can dominate in the enhanced NS cooling as well provided the nucleons of one species 
are strongly superfluid while the other ones are moderately superfluid. Very strong superfluidities of 
neutrons and protons (upper right corners of the figures) switch off all neutrino processes involving 
nucleons. As a result, the neutrino bremsstrahlung due to electron-electron collisions, which is prac- 
tically unaffected by superfluidity, becomes dominant. It is difficult to expect that this weak process 
is important under other conditions. 
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7 Cooling of neutron stars 

7.1 Review of articles on neutron star cooling 

7.1.1 Overall review 

The theory of NS cooling has been developing over more than 30 years. The first articles appeared even 
before the discovery of NSs. Their authors tried to prove that not too old NSs may emit sufficiently 
powerful thermal X-ray radiation which could serve to discover NSs. The first estimates of the thermal 
emission from cooling NSs were most probably done by R. Stabler (1960). Four years later Chiu (1964) 
repeated these estimates and theoretically proved the possibility to discover NSs from their thermal 
emission. First, simplified calculations of the NS cooling were done by Morton (1964), Chiu and 
Salpeter (1964) and also by Bahcall and Wolf (1965b) after the discovery of X-ray sources in the Crab 
nebula and Scorpion constellation in the balloon experiments by Bowyer et al. (1964). 

The foundation of the strict cooling theory was laid in the fundamental paper by Tsuruta and 
Cameron (1966) who explicitly formulated the main elements of the theory — the relationship between 
the internal and surface NS temperatures, the neutrino and photon cooling stages, etc. Later the theory 
has been developed in many papers. 

In 1970s and 1980s, the main attention was paid to the equation of state and nuclear composition 
of the NS cores (including possible appearance of exotic particles), to neutrino reactions and to the 
relationship between the internal and surface temperatures of NSs (particularly in the presence of 
strong magnetic fields). The achievements of the theory in the middle of 1970s were reviewed by 
Tsuruta (1979). A splash of theoretical activity in the beginning of 1980s was concerned with the 
launch of the Einstein space observatory (Sect. 8.1). The articles of the "Einstein series" were also 
reviewed by Tsuruta (1980, 1981, 1986). The launch of the ROSAT observatory in 1990 (Sect. 8.1) 
initiated a new rise of theoretical activity which is being continued up to now. Recent articles are 
mainly focused on the mechanisms of NS reheating at the late evolutionary stages (see below) and 
the effects of superfiuidity in the NS cores. Let us mention recent review by Tsuruta (1998). New 
microscopic theories of dense matter are appearing permanently complicating cooling models; many 
problems are still unsolved. 

Let us outline some aspects of the theory. 

The first articles were devoted to the standard cooling (produced by the neutrino reactions from 
the "standard" collection. Sect. 5). The enhanced cooling has been simulated since the end of 1970s. 
It has been assumed in many articles that the enhanced neutrino luminosity is associated with exotic 
composition of the NS cores containing pion condensates of quark plasma; see, e.g., Tsuruta (1979), 
Glen and Sutherland (1980), Van Riper and Lamb (1981), Yakovlev and Urpin (1981), Richardson 
et al. (1982), Umeda et al. (1993, 1994a, b), Schaab et al. (1996). By the end of 1980s a new exotic 
cooling agent, kaon condensate, has been introduced (Page and Baron, 1990; Schaab et al., 1996; Page, 
1997). New models of dense matter with highly polarized pion degrees of freedom (Voskresensky and 
Senatorov, 1984, 1986; Schaab et al., 1997b) have been proposed. According to these models the 
neutrino luminosity is strongly enhanced due to virtual excitation of pion field even if density of 
matter does not exceed the critical density of actual pion condensation. Modern cooling theories of 
stars with the quark cores, pion or kaon condensates have been described, for instance, by Pethick 
(1992) and Schaab et al. (1996). 

A new stage of the theory of enhanced cooling has been opened by Lattimer et al. (1991) who 
have shown that the direct Urea process can be allowed in the NS cores with the standard nuclear 
composition for many realistic equations of state. The process initiates rapid NS cooling (e.g.. Page 
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and Applegate, 1992; Van Riper and Lattinier, 1993; Gnedin and Yakovlev, 1993; Page, 1994; Gnedin 
et al., 1994; Page, 1995b; Levenfish and Yakovlev, 1996; Schaab et al., 1996) without invoking "ex- 
otic" hypotheses (Sect. 2.2). Detailed description of the "non-exotic" standard and enhanced cooling 
theories has been given by Pethick (1992). 

The effect of the NS magnetic field on the relationship between the internal and surface stellar 
temperatures has been taken into account starting from the articles by Tsuruta and coauthors (Tsu- 
ruta et al., 1972; Tsuruta 1974, 1975). These articles were the first where a representative set of 
neutrino reactions in the NS crust and core were included (plasmon decay, annihilation of electron- 
positron pairs, photon decay, electron bremsstrahlung due to scattering off nuclei, neutron branch of 
the modified Urea process). A detailed study of the relationship between the internal and surface 
temperatures in a NS with the magnetic field normal to the surface was carried out by Van Riper 
(1988). He also analysed in detail (Van Riper, 1991a) the effect of such magnetic fields on the NS 
cooling. Page (1995a) as well as Shibanov and Yakovlev (1996) considered the cooling of a NS with 
a dipole magnetic field and showed that the dipole field affected the cooling much weaker, and in a 
qualitatively different way, than the purely radial magnetic field. 

Let us mention also a recent series of articles by Heyl and coauthors (Heyl and Hernquist 1997a, b, 
1998a, b, c; Heyl and Kulkarni, 1998) devoted to the cooling of NSs with superstrong magnetic fields 
IQ-'^^-lO^^ G (the so called "magnetars", see Sect. 8.1.1). These fields may strongly reduce thermal 
insulation of the NS envelope, making the magnetar's surface much hotter at the early cooling stage, 
than the surface of an "ordinary" NS. Let us stress that microscopic properties of matter in superstrong 
magnetic fields (equation of state, thermal conductivity) are poorly known, so that the results by Heyl 
and coauthors can be regarded as very preliminary. 

The cooling can also be affected noticeably (Chabrier et al., 1997; Potekhin et al., 1997; Page, 
1997, 1998a, 1998b) by the presence of a thin (of mass < W~^Mq) envelope of light elements (H, He) 
at the surface of a non-magnetized or weakly magnetized NS. Owing to the higher electron thermal 
conductivity of plasma composed of light elements, the NS surface appears to be significantly warmer 
at the early cooling stages. 

An important contribution into the theory was made in the PhD thesis by Malone (1974). He 
was the first who calculated the standard NS cooling beyond the approximation of isothermal internal 
layers. This allowed him to describe thermal relaxation of the internal layers in the first 100-1000 
years of NS life. Analogous nonisothermal calculations of the cooling enhanced by the presence of the 
pion condensate were done in the PhD thesis of Richardson (1980). The results of both theses were 
published in one article (Richardson et al., 1982). Later the thermal relaxation of the internal layers 
in a young NS was studied in a series of articles by Nomoto and Tsuruta (1981, 1982, 1983, 1986, 
1987) as well as by Lattimer et al. (1994). The thermal relaxation is accompanied by propagation of a 
cooling wave from the internal layers to the surface. In principle, the appearance of the cooling wave 
at the surface can be observed in young NSs. The moment of appearance depends on the equation of 
state in the central stellar region (Lattimer et al., 1994). 

In 1980 Glen and Sutherland and a year later Van Riper and Lamb (1981) included the effects 
of General Relativity into the equations of NS thermal evolution (in addition to the equations of 
hydrostatic equilibrium as had been done before) . In this connection let us mention a recent article by 
Schaab and Weigel (1998) who carried out the first two-dimensional cooling calculations of a rotating 
NS with exact account of the general relativistic effects produced by rotation. 

New direction of study was opened by Alpar et al. (1987) and Shibazaki and Lamb (1989). The 
authors took into account possible reheating at late cooling stages (NS age t > 10^ yr) due to viscous 
dissipation of rotational energy inside a NS. The effect is caused by interaction of superfluid and 
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normal components of matter in the inner crust of a pulsar which is spinning down under the action 
of magnetodipole losses. The cooling theory with viscous reheating has been developed further in a 
number of articles (e.g., Van Riper, 1991b; Van Riper et al., 1991, 1995; Umeda et al., 1993, 1994b; 
Shibazaki and Mochizuki, 1995; Page, 1997, 1998a, 1998b). 

The core of a cooling NS can also be reheated by the energy release associated with a weak deviation 
from beta-equilibrium (Reisenegger, 1995). In addition, the reheating of the star with a non-superfluid 
core can be produced by ohmic dissipation of the core magnetic field (Haensel et al., 1990; Yakovlev, 
1993; Shalybkov, 1994; Urpin and Shalybkov, 1995) due to the enhancement of the electric resistance 
across the strong magnetic field. Criticism of this effect by Goldreich and Reisenegger (1992) is not 
convincing because, while calculating the electric resistance, the latter authors neglected motion of 
neutron component of matter (inconsistency of such approximation is clearly seen from the results 
by Yakovlev and Shalybkov, 1991). Finally, the reheating in an old NS (t > lO''' yr) can be provided 
by ohmic decay of the magnetic field in the NS crust (Miralles et al., 1998; Urpin and Konenkov, 
1998). Let us add that the magnetic field dissipation in the crust represents a very important process 
which determines evolution of the surface magnetic field, magnetodipole NS spindown, NS activity 
as radio pulsar, etc. The dissipation process is associated with cooling since the electric resistance of 
the curst decreases in the course of cooling (until it reaches minimum produced by electron scattering 
off charged impurities). Therefore, the dissipation rate depends on cooling type: it is very weak for 
the rapid cooling. On the other hand, the magnetic field decay weakens magnetodipole spindown 
of the star. This means that cooling, magnetic field evolution and spindown should be considered, 
in principle, selfconsistently as magneto-rotational evolution of a NS. Corresponding theory has been 
developed in a series of papers by Urpin and coauthors (see, e.g., Urpin and Muslimov, 1992; Urpin 
and Van Riper, 1993; Urpin and Konenkov, 1998, and references therein) without account for the 
effect of superfiuidity on the NS cooling. 

7.1.2 Cooling of superfluid stars 

The effect of nucleon superfluidity in the NS core on the cooling has been taken into account starting 
from the article by Tsuruta et al. (1972). In the first articles, as a rule, it has been included in a sim- 
plified manner; it has been assumed that the superfiuidity completely switches off corresponding heat 
capacity or neutrino emission from the very beginning of the cooling (which corresponds to infinitely 
large critical temperatures). Maxwell (1979) was the first who studied the dependence of the standard 
NS cooling on nucleon critical temperatures. Friman and Maxwell (1979) considered a simplified NS 
model with constant-density core neglecting the effects of General Relativity, but implying sufficiently 
realistic neutrino luminosity. For a long time, before publication of the article by Page and Applegate 
(1992), the superfluidity was not considered as a powerful cooling regulator. It was taken into account 
because it appeared inevitably in microscopic theories of dense matter (Sect. 3.1). However recent 
studies show that the superfluidity is one of the most important factors which affect the standard and 
enhanced NS cooling. 

Simulations of the standard cooling of NSs with superfluid cores have been performed by 
many authors, in particular, by Tsuruta et al. (1972), Malone (1974), Tsuruta (1978, 1979, 1980), 
Maxwell (1979), Glen and Sutherland (1980), Van Riper and Lamb (1981), Nomoto and Tsuruta 
(1981, 1986, 1987), Richardson et al. (1982), Page (1994), Levenfish and Yakovlev (1996), Schaab et 
al. (1996), Levenfish et al. (1999), Yakovlev et al. (1999). While calculating the standard neutrino 
luminosity of the NS core, as a rule, one took into account the contribution of the neutron branch of 
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Figure 14: Cooling curves calculated using accurate (solid and dash-dotted lines) and approximate (dashed and dotted 
lines) description of the effects of superfluidity on the neutrino luminosity and heat capacity (see text). Solid and dashed 
lines correspond to the enhanced cooling of a NS (M = 1.48 M©) in the presence of the proton superfluidity with 
Tcp = 10* K. Dash-dotted and dotted lines refer to the standard cooling of a NS (1.30 Mq) in the presence of neutron 
and proton superfluidities with Ten — 10^ K and Tcp = 10* K. The parameters of the models are given in Sect. 7.2.2. 
Neutrino emission due to nucleon pairing is neglected. 

the modified Urea processes and the neutrino bremsstrahlung due to nn- and np-scattering. Tsuruta 
et al. (1972) were the first who studied the NS coohng accompanied by the singlet-state pairing of 
neutrons and protons. The authors described reduction of the n-branch of the modified Urea-process 
and the wp bremsstrahlung by the asymptotic expressions obtained by Itoh and Tsuneto (1972) for 
the modified Urea in the limit of strong (T <C T,^ singlet-state pairing of neutrons; they multiplied the 
nn-scattering rate by exp(— A„/A;bT). Maxwell (1979) just switched off the modified Urea process at 
temperatures below T^n or Tcp. In a number of articles various authors used the simplified reduction 
factors proposed by Malone (1974): exp[— (A^ -|- Ap)/A;bT] for superfluidity of neutrons and proton 
at once; exp(— Atv/^bT) for singlet-state superfluidity of neutrons or protons (for the neutron branch 
of the modified Urea-process and np-bremsstrahlung) . The reduction of the nn-bremsstrahlung by 
the neutron superfluidity was described by the expression exp(— 2A„/A:bT'). The accuracy of these 



approximations was analysed in Sects. 4 and 5. Their effect on the cooling is demonstrated in Fig. 14 



The complete set of neutrino reactions considered in Sects. 4 and 5 was first used by Levenfish and 
Yakovlev (1996). These authors implied also a more accurate description of the effects of superfluidity 
on the heat capacity and neutrino luminosity (Sects. 3-5), but neglected the neutrino emission due 
to Cooper pairing of nucleons (Sect. 6). The latter process has been included only in the recent 
simulations by Schaab et al. (1997b), Page (1998a, 1998b), Levenfish et al. (1998, 1999), and Yakovlev 
et al. (1998, 1999). Let us mention that Levenfish et al. (1998, 1999) and Yakovlev et al. (1998, 1999) 
have made use of more accurate expressions for the neutrino emissivity given in Sect. 6 (although 
neglected the relativistic correction to the neutrino emissivity due to pairing of protons). 
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Figure 15: Standard cooling of a NS with M = 1.30 M© (Fig. a) and enhanced cooling of a NS with M = 1.48 M© 
(Fig. b) in the presence of superfluidity. The curves are labeled by the values of Ig Ten and Ig Tcp. Solid lines are obtained 
including the neutrino emission due to pairing of nucleons while dashed lines are obtained neglecting this emission. 
Dotted lines show cooling of non-superfluid NSs. The parameters of the models are given in Sect. 7.2.2. The solid and 
dashed lines for IgTcn = 7.8 and IgTcp = 7.8 in Fig. (a) coincide with each other and so do the lines for IgTcn = 8.3, 
IgTcp = 8.6 and IgTcn = 10.0, IgTcp = 8.4 in Fig. (b). 

The importance of neutrino emission due to Cooper pairing of nucleons for the standard NS coohng 



is illustrated in Fig. 15 a. One can see that the appearance of superfluidity at the neutrino cooling stage 
can not only slow down the NS cooling, as thought before, but, on the contrary, strongly accelerate it 
(see, e.g., Yakovlev et al., 1999). The acceleration can be so fast that at t ~ 10^-10^ yr the star can be 
much cooler that in the case of the enhanced neutrino energy losses. An exclusion is provided by super- 
fluidity of protons alone since proton pairing does not produce intense neutrino emission. In addition, 
the difference between the cooling curves is small in the case of normal protons and strongly superfluid 
neutrons {Ten ^ 10^ K), in which the Cooper pairing of neutrons is not a dominant process (Sect. 6.2). 



Simulation of enhanced cooling of NSs with the standard composition of the cores was started 
by Page and Applegate (1992). These authors were the first who discovered an interesting feature in 
the cooling of NSs in which the direct Urea process was open and neutrons or protons were superfluid: 
after thermal relaxation was over the surface temperature of these stars fell down rapidly to the value 
Tg = Ts{Ti), appropriate to the internal temperature Tj = aTc, and remained almost constant during 
subsequent neutrino cooling stage. Here Tc is the critical temperature of neutrons or protons in the 
core, and a ~ 1 is a numeric coefficient (for instance. Page and Applegate, 1992, obtained a ~ 0.2). 
Moreover, the temperature Tj is almost insensitive to other NS parameters (Page, 1994). The nature 
of this phenomenon is very simple — the star cools rapidly before the superfluidity onset, whereas 
the superfluidity suppresses the cooling and "freezes" the internal temperature at the level of aTc. 
This opens possibility to study superfluidity in those NSs whose surface temperature is known from 
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observations. 

Page (1994) analyzed this possibility in more detail (neglecting the Cooper-pairing neutrino emis- 
sion). He showed that the measured surface temperatures of PSR 0656+14 (Finley et al., 1992) and 
Geminga (Halpern and Ruderman, 1993) could be explained in the standard and enhanced cooling 
models. The observations of the Vela pulsar (Ogelman et al., 1993) are more easily explained by 
the enhanced cooling model, and the observations of PSR 1055-52 (Brinkmann and Ogelman, 1987) 
by the standard cooling. In all these cases one needs (Page, 1994) the presence of the neutron and 
proton superfluidities with high critical temperatures Tc ~ 10^ K in the entire NS body where the 
main neutrino emission occurs. 

Let us stress that the results and quantitative conclusions by Page and Applegate (1992) and 
Page (1994) are considerably modified by including the effects of joint superfluidity of nucleons and 
Cooper-pairing neutrinos. However the principal conclusion by Page and Applegate remains the same: 
a NS of age 10^-10^ yr with superfluid nucleons in its core is a "thermometer" of this superfluidity 
— one can measure (constrain) the critical temperatures of neutrons Ten and protons Tcp from the 
values of the surface temperature. 

Aside of Page and Applegate (1992) and Page (1994) , the enhanced cooling of NSs with the 
standard nuclear composition and nucleon superfluidity in the NS cores has been studied, for instance, 
by Gnedin and Yakovlev (1993), Van Riper and Lattimer (1993), Gnedin et al. (1994), Levenfish and 
Yakovlev (1996), Page (1998a, b), Levenfish et al. (1998, 1999), Yakovlev et al. (1998, 1999). As 
in the case of the standard cooling, the effect of superfluidity on the heat capacity and neutrino 
luminosity was described in earlier articles by the approximate factors of the form ex.p{— A^ / kBT) , 
but gradually the theory presented in Sects. 3-5 has become implemented. The difference of the 
standard and enhanced cooling curves calculated using the accurate (Sects. 3-5) and approximate 
(exp(— Aj//cbT)) reduction factors is illustrated in Fig. 14. The accurate reduction factors "freeze" 



the cooling at the neutrino cooling stage at a lower internal temperature T and, accordingly, at a 
lower surface temperature (Gnedin and Yakovlev, 1993). 

In calculations of the enhanced neutrino luminosity the contribution from Cooper pairing of nucle- 
ons (Sect. 6) has usually been neglected (excluding the articles by Schaab et al., 1997b; Page 1998a, 
b; Levenfish et al., 1998, 1999; Yakovlev et al., 1998, 1999). According to Yakovlev et al. (1999) 
this process is especially important in the presence of the proton superfluidity with Tcp ^ Tc„. Such 
superfluidity appears at the early cooling stage and suppresses the direct Urea process long before the 
onset of the neutron pairing. Therefore, a splash of the neutrino emission associated with the neutron 
pairing is very pronounced. In this case the NS cooling does not slow down, as it would be without 
neutrino emission due to nucleon pairing, but is strongly accelerated (Fig. p^b). Let us remind that 
similar situation takes place in the case of the standard cooling (also see Sect. 7.2). 

As a result, the nucleon superfluidity plays the leading role for the both, the standard 
and enhanced, cooling (Page 1998a, b; Levenfish et al., 1998, 1999; Yakovlev et al., 1998, 1999). 
Moreover, the superfluidity mixes the cooling types, since the enhanced cooling may look like the 
standard one while the standard cooling may look like enhanced (especially if the neutrino emission due 



to Cooper pairing is taken into account, cf. Figs. |l5|a and |l5p). In these cases, the star, independently of 



its cooling type, remains to be a good "thermometer" for measuring the nucleon critical temperatures: 
even weak variations of Ten or Tcp affect strongly the cooling curves. The only problem is to "calibrate" 
properly this thermometer (Sects. 3-6). 

Finally, let us mention that Schaab et al. (1998b) have considered recently, for the first time, the 
standard and enhanced cooling in the presence of the triplet-state neutron superfluidity of type C 
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(with nodes of the gap at the Fermi surface). The authors have been based on the model calculations 
of neutron pairing in the strong magnetic fields by Muzikar et al. (1980) according to which the field 
B > 10^^ G makes the pairing of type C energetically more favourable than that of type B. Let us 
notice also the first calculations of cooling (Schaab et al., 1998a) of a NS with the hyperonic core taken 
into account the hyperon superfluidity. We mention, however, that Schaab et al. (1998a, 1998b) have 
made use of approximate reduction factors of neutrino reactions and approximate expressions for the 
neutrino emissivity due to Cooper pairing of particles. 

The theoretical conclusion on the leading role of superfluidity in NS cooling was made at about 
the same time when new observational data on the surface temperatures of some NSs appeared (Sect. 
8) . The latter data brought the problem of studying superfluidity in the NS cores to reality. 

7.2 Cooling simulations 

In this section we illustrate the possibility to explore nucleon superfluidity in the NS cores. We will 
mainly follow our recent article (Levenfish et al., 1999). In addition, we include into the simulations 
the relativistic correction to the neutrino luminosity produced by proton pairing (Sect. 6.1). It affects 
slightly the cooling curves at some values of Ten and Tcp but does not change the main conclusions by 
Levenfish et al. (1999). 

7.2.1 Simulation scheme 

Below we present the results obtained with the cooling code (Gnedin and Yakovlev, 1993; Gnedin 
et al., 1994; Haensel and Gnedin, 1994; Levenfish and Yakovlev, 1996; Levenfish et al., 1998) that 
is based on the approximation of isothermal stellar interiors. The code was mainly constructed by 
O.Y. Gnedin. The isothermal approximation is valid for a star of age t > (10-10^) yr, inside which 
the thermal relaxation is over. The superfluidity affects strongly the surface temperature only after 
reaching the thermally relaxed stage. Following Glen and Sutherland (1980) we have assumed that 
the isothermal region is restricted by the condition p > p5 = lO^'^ g cm~^. The real boundary of the 
isothermal region is situated at lower p and depends on temperature: it moves to the surface in the 
course of cooling. The chosen value of ph guarantees that the region p > pi, is isothermal in a star of 
age t > 10'^ yr. The quantity Tj(t) = T{r,t) exp[<l>(r)] (which may be called the internal redshifted 
temperature) is constant over the isothermal region at any moment of time t. Here T{r, t) is the local 
temperature, $(r) is the dimensionless gravitational potential, and r is radial coordinate. 

We take into account explicitly the effects of General Relativity on the NS structure and cooling. 
Cooling simulation in the isothermal approximation is reduced to solving the equation of thermal 
balance (see, e.g.. Glen and Sutherland, 1980) 

An~' 
cm) -^ = -L,i{Ti) - L^,{Ts) ; (85) 

// 2Cm \ ^^/^ 
c^dy, dy = 47rr2 f 1 \ dr , (86) 

Lui= I Q exp(2$) dy, L^i = 47r7^V^/ exp(2$b). (87) 

In this case, C is the total heat capacity of the star, c^, is the specific heat capacity, and ra = m[r) is 
the gravitational mass inside a sphere with radial coordinate r. The quantities Lj^i and L^i determine 
the neutrino and photon NS luminosities, respectively, Q is the neutrino emissivity, $5 is the value 
of ^ at the isothermal region boundary {p = ph), and a is the Stefan-Boltzmann constant. 
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The photon luminosity of the star depends on the effective temperature Tg of its surface. The 
relationship between Tg and the temperature Tj = Tjexp(— <I>t) at p = Pb is determined by thermal 
insulation of the outer envelope {p<pb), where the main temperature gradient is formed. We have not 
taken into account the effect of the magnetic field on the NS cooling and have used the relationship 
Ts = Ts{Tb), obtained recently by Potekhin et al. (1997) for B = 0. We assume that the NS surface 
may be covered by a thin layer (of mass < 10^^^ -^o) of hydrogen or helium. This amount of light 
elements is too small to affect the thermal insulation of the envelope and the NS cooling but it can 
affect the spectrum of thermal radiation. Actually, the effect of the dipole surface magnetic fields 
B < 5 X 10^^ G on the NS cooling is rather weak (Shibanov and Yakovlev, 1996). Therefore, our 
calculations can be used, at least qualitatively, for the stars with such fields. In these cases, under Tg 
we mean the average effective surface temperature which determines the total (nonredshifted) photon 
surface luminosity of the star as L^ = AiraTZ'^T^. 

7.2.2 Models of cooling neutron stars 

For simplicity, we assume that matter of the NS core consists of neutrons, protons and electrons (muons 
and hyperons are absent). We will adopt a moderately stiff equation of state proposed by Prakash 
et al. (1988) (the version with the compression modulus Kq = 180 MeV and the same simplified 
symmetry factor Sy, that was used by Page and Applegate, 1992). The maximum NS mass, for a 
given equation of state, is 1.73 Mq. In order to study the enhanced and standard cooling we consider 
the NS models of two masses. In the first case, the NS mass is M = 1.48 Mq, radius TZ = 11.44 km, 
and the central density /Jc = 1.376x 10^^ g cm~^, while in the second case M = 1.30 Mq, TZ = 11.87 km, 
and pc = 1.07 x 10^^ g cm~^. The adopted equation of state opens the direct Urea process at the density 
/J > Per = 1-30 X 10^^ g cm~^. Therefore, the 1.48 Mq NS suffers the enhanced cooling: the powerful 
direct Urea process (Q) is allowed in a small central kernel of radius 2.32 km and mass 0.035 Mq (in 
addition to the processes (|5)"(^ and (^) in the entire core). In the 1.30 Mq NS the threshold value 
Per is not reached, and the star has the standard neutrino luminosity, determined by the processes 
(jl|)"(^) and (|5|). Notice that Levenfish and Yakovlev (1996), while calculating the equation of state, 
set the parameter no (the standard saturation nuclear-matter density) equal to 0.165 fm^'^. In the 
present calculations we set no = 0.16 fm~'^. Owing to this reason the mass of the rapidly cooling NS 
model is somewhat different from that (1.44 Mq) used by Levenfish and Yakovlev (1996). 

In order to calculate the heat capacity of the superfluid NS core and its neutrino luminosity due to 
the reactions (|l])-(|5|) we have used equations of Sects. 3-6. The results obtained in these sections are 
summarized in Sect. 6.2. The neutrino bremsstrahlung due to ee scattering in the NS core has been 
neglected: we have shown that its effect on the cooling is negligible, for our NS models. We have also 
taken into account the neutrino luminosity of the NS crust due to bremsstrahlung of electrons which 
scatter off atomic nuclei (using an approximate formula proposed by Maxwell, 1979): Lbr = 1.65 x 10'^^ 
{Mcr/MQ){Tb/W^Kfex.p{2^b) erg s-\ where Mcr is the crust mass. In our case, M^r = 0.0120 Mq 
for the 1.48Mq model, and Mcr = 0.0153 Mq for the 1.3 Mq model. The NS heat capacity has been 
set equal to the sum of the partial heat capacities of n, p and e in the stellar core; the contribution 
of the crustal heat capacity has been neglected due to the smallness of the crust mass for the chosen 
NS models. In the calculations of the neutrino luminosity and the heat capacity the effective masses 
of neutrons and protons in the NS cores have been taken as m'^ = 0.7 mjv 

Various microscopic theories of Cooper pairing in the NS cores give very wide scatter of the critical 
temperatures of the neutron and proton superfluids , Ten and Tcp (~ 10^-10^" K), and various density 
dependence of these temperatures (Sect. 3.1). Thus we have made a simplified assumption that the 
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critical temperatures T^n and T^p are constant throughout the NS core and may be treated as free 
parameters. This is the main simphfication of our simulations. We assume that protons suffer ^So-state 
pairing (superfluidity of type A, Table 1) and neutrons suffer ^P2-state pairing with zero projection 
of the pair moment onto quantization axis (superfluidity of type B). We will study the cooling for 
different values of T^n and T^p, and we will determine (Sect. 8.2.2) those values which are in better 
agreement with observations. 

As already mentioned in Sect. 7.1, the difference between the enhanced and standard cooling 
regimes is smeared out in the presence of superfluidity. Under these conditions, the cooling regime 
formally indicates if the direct Urea is allowed or forbidden in the NS core. 

7.2.3 Results 

We have calculated about 2000 NS cooling curves. The curves describe the dependence of the effective 



surface stellar temperature T^ = TgJl — TZg/TZ, as detected by a distant observer, on age t; here, 
Ug is the gravitational NS radius. For M = 1.30 Mq and M = 1.48 M© we have T^ /Ts = 0.822 and 
0.786, respectively. The critical temperatures of neutrons and protons, T^n and T^p, in the NS core 
have been varied in a wide interval from 10^ to lO^'^ K. 

The examples of the cooling curves for selected values of T^n and Tf.p were already discussed above 
(Figs, [iq a and |l^b). However it is rather inconvenient to analyse the results in the form of the cooling 
curves. It is better to plot the values of (Ten, Tcp), which lead to certain surface temperatures T"^ of 



the NS of given age t. Figures 16- 1£ present the results in this way. For example, we have chosen the 
values of t, which correspond to the ages of the Geminga pulsar (3.4 x 10^ yr. Figs. 16 and [T^ and 



PSR0656-I- 14 (10^ yr. Figs. 17 andH). Observational data on these and other cooling NSs are given 
in Sect. 8.1. 

Before discussing the results of calculations let us remind some relationships which determine the 



main features of the NS cooling. We will use these relationships for describing Figs. 16-19- 

The effects of the neutron and proton superfluidities on the heat capacity are different. In a non- 
superfluid NS, the heat capacities oip and n constitute ~ 1/4 and ~ 3/4 of the total heat capacity Ctot^ 
respectively. Therefore, a strong superfluidity of p reduces Ctot by ~ 25%, while a strong superfluidity 
of n reduces it by about 4 times. When the ratio of the internal stellar temperature Tj to a critical 
nucleon temperature Tcn decreases from 0.3 to 0.1 (for instance, in the course of NS cooling), the 
heat capacity of nucleons N = n oi p decreases by more than three orders of magnitude and becomes 
much lower than the heat capacity of electrons. Further reduction of the nucleon heat capacity does 
not change the total heat capacity. The appearance of the weak superfluidity of nucleons A^ almost 
doubles their heat capacity due to the latent heat release at the phase transition. The heat capacity of 
nucleons remains higher than the heat capacity of non-superfluid nucleons as long as Ti/TcN decreases 
from 1 to 0.3, i.e., as long as the difference of temperature logarithms is \gTcN — Ig^j < 0.5. 

The effects of superfluidities of n and p on the NS neutrino luminosity are also different. For the 
direct and modified Urea processes this difference is weak. If Cooper-pairing neutrino emission were 
neglected, the asymmetry of Figs. p!^19 with respect to inversion of the axes T^n ^ ?cp, would be 



mainly explained by different contributions of n and p into the heat capacity. This can be proved by 



comparing Figs. 16-19| with analogous figures in Levenfish and Yakovlev (1996), where the neutrino 
emission due to Cooper pairing was neglected. Inclusion of the Cooper-pairing neutrinos greatly am- 
plifies the asymmetry: the neutrino emission due to proton pairing is weak, while the emission due 
to neutron pairing dominates over the standard neutrino losses for T, < 10^ K < T^n and over the 
direct Urea process for Ti < T^n 'C T^p- The main neutrino energy release in the process takes place 
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Figure 16: Lines of values of Ten and T^p which correspond to certain internal temperatures Ti (the values of IgTi are 
given near the curves) or surface temperatures T^ (given in Table 0) of a NS with the enhanced neutrino luminosity 
(1.48 M0) and the Geminga's age (3.4 ■ 10^ yr). The darker the background the cooler is the star. The region of joint 
neutron and proton superfluidity (in the center and the upper right corner) is enclosed by the dashed lines. A small 
region where the superfluidity does not appear by the given age (the left lower corner) is separated by the thick solid 
line. 

at 0.2 < Ti/Tcn < 0.96 (as long as the difference of the temperature logarithms is IgT^ — Ig^i < 0.7). 

Enhanced cooling. Figures O and O illustrate the enhanced cooling of NSs of the ages of the 
Geminga and PSR0656 + 14 pulsars. 



Figure 16 shows isotherms of the internal temperature Tj of the NS of the Geminga's age versus 



Ten and Tcp. Since the temperature is related to the surface temperature, the same lines are isotherms 
of the surface temperature T^ (Table |2|). Dashes exhibit the auxiliary lines which enclose the region 
of the joint superfluidity of nucleons. The region to the left of the upper dashed line corresponds to the 
superfluidity of protons and normal neutrons, while the region to the right of the lower line corresponds 
to the superfluidity of neutrons and normal protons. The dashed lines intersect at the isotherm of 
the temperature Tj = 10^'^ K, to which a non-superfluid star would cool down by the moment t. 
This isotherm, plotted by the thick line, encloses the region where the nucleon superfluidity has not 
appeared by the moment t and does not affect the cooling. Notice, that owing to the effects of General 
Relativity, the isotherms Tj correspond to somewhat higher local temperatures of matter (see above). 
Thus, the isotherm Tj = 10^'^ K is associated with the values of T^n which are slightly higher than Tj. 
First we discuss behaviour of isotherms to the right of the lower auxiliary line. The vertical 
segments of the isotherms reveal that the NS cooling is governed by the only superfluidity of n. With 
increasing Ten at Ig Tc„ > 7.7, the core temperature Tj grows up which can be explained like this. 
The higher T^n, the earlier the neutron superfluidity appears. Accordingly, the powerful direct Urea 
process is suppressed earlier, cooling delay is longer, and the NS is hotter at given t. For low Ten, the 
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Figure 17: Same as in Fig. ^ but for a NS of the PSR0656 + 14 age (10^ yr). 

NS cools down in a non-super fluid state almost all its life and is sufficiently cold by the age t. The 
neutrino luminosity of such a star is rather weak and becomes comparable to the photon luminosity. 
Additional neutrino energy losses due to neutron pairing at IgTc^ > 6.6 slightly accelerate the cooling 
at 6.6 < IgTcn < 6.6 + 0.7. However, for IgTcn ^ 6.6 + 0.5 the acceleration is compensated by the 
latent heat release produced by the jump of the neutron heat capacity. The values 0.7 and 0.5 have 
been explained above. 

Vertical segments of the isotherms intersect the auxiliary line at temperatures Tcp, at which the 
proton superfluidity appears. With further increase of Tcp, the NS heat capacity (determined by 
protons for IgTcn > 7.4) jumps up by about twice. The latent heat release increases Tji in a strip of 
width ~0.5 in IgTcp above the lower dashed line the isotherms shift to the left. With further growth 
of Tcp the heat capacity is strongly reduced, the NS becomes cooler and the isotherms shift to the 
right. 

In a similar fashion, exchanging the proton and neutron superfluids, we can explain horizontal 
segments of the isotherms to the left of the upper auxiliary line. In this case acceleration of the NS 
cooling in the range IgTcp > 6.6 does not appear because of the weakness of the neutrino emission 
produced by the proton pairing. A strong proton superfluidity suppresses the heat capacity weaker 
than a neutron superfluidity with the same critical temperature. Therefore, at high Tcp the star is 
warmer than for the same Ten at vertical segments of isotherms in the lower right part of the figure. 

The horizontal segments of the isotherms intersect the upper auxiliary line at temperatures Ten 
which switch on the neutron superfluidity. This superfluidity induces the latent heat release and 
associated slight NS heating (dips on the isotherms to the right of the auxiliary line). With increasing 
Tcp, the NS heating produced by the neutron superfluidity becomes weaker and disappears (the dips 
vanish). The effect is mainly caused by the neutrino emission due to neutron pairing which is more 
pronounced if the direct Urea process is strongly suppressed. At high Tcp the neutron pairing becomes 
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Table 2: Relationship between Tj and T^ [K] for Fig. |l§^ 



M = 1.48 Mq 


igr, 
igrr 


5.10 5.30 5.50 5.70 5.90 6.05 6.20 6.36 6.40 6.44 6.57 6.60 6.61 6.65 6.73 
4.53 4.61 4.68 4.76 4.84 4.91 4.98 5.05 5.07 5.09 5.16 5.17 5.18 5.20 5.24 


igr, 
igrr 


6.80 6.89 7.00 7.07 7.20 7.25 7.40 7.43 7.60 7.61 7.77 7.79 7.87 7.97 8.10 
5.28 5.32 5.38 5.42 5.49 5.52 5.60 5.62 5.71 5.72 5.81 5.82 5.86 5.92 5.99 



M = 1.3OM0 


Igr, 

igrr 


5.73 5.75 5.83 5.90 5.98 6.30 7.00 7.21 7.23 7.29 7.32 7.34 
4.78 4.79 4.82 4.85 4.88 5.02 5.37 5.48 5.49 5.53 5.54 5.55 


Igr, 

ig^r 


7.38 7.43 7.52 7.54 7.58 7.65 7.67 7.70 7.75 7.79 7.80 7.81 
5.57 5.60 5.65 5.66 5.69 5.72 5.74 5.75 5.78 5.80 5.81 5.81 


Igr. 

igrr 


7.83 7.85 7.86 7.88 7.89 7.98 7.99 8.01 8.05 8.06 8.12 8.14 
5.82 5.84 5.84 5.85 5.86 5.91 5.91 5.93 5.95 5.95 5.99 6.00 



the main neutrino emission mechanism. Since the neutrino emission due to the neutron pairing is 
more efficient than the latent heat release the appearance of the neutron superfluidity with the growth 
of Ten does not delay the cooling but, on the contrary, accelerates it. A noticeable NS cooling via 
Cooper-pairing neutrinos takes place in a strip of the width 0.7 (in IgTcn) to the right of the upper 
dashed line. Lower and to the right of this strip the cooling is associated with reduction of the neutron 
heat capacity. 

Now consider the region of Ten, T^p between the auxiliary lines in more detail. Some increase of Tj 
at Ig Ten ^ 6.6 + 0.5 and Ig T^p < 6.6 + 0.5 is caused by the latent heat release at onset of superfluidity 
of n and p. Further growth of Ten ^ T^p induces initially rapid decrease and then weak increase of Tj. 
The decrease is explained by exponential reduction of the heat capacity by the joint superfluidity of 
n and p while the weak increase is associated with reduction of the direct Urea process at the early 
cooling stages. 

Figure ^ is analogous to Fig. |l^, but corresponds to a younger NS of age t = 10^ yr. Figures |l^ 

or Ten *C Tcp to support high surface temperature Tf° (or Tj) 



and 17 show that one needs T^p <^ T^ 



for a longer time. 



Standard cooling. Standard cooling of a 1.30 Mq NS is illustrated in Figs. |18| and |19| for the stars 
of Geminga and PSR 0656 + 14 ages, respectively. Isotherms are qualitatively different from those for 



the enhanced cooling (cf. with Figs. |l^ and 17): even an approximate symmetry of the neutron and 
proton superfluidities is absent. The asymmetry can be attributed to the weakness of the standard 
neutrino energy losses. First, in the absence of such a powerful cooling regulator as the direct Urea 
process the difference of the heat capacities of n and p (see above) is more pronounced. Second, 
the Cooper-pairing neutrino emission becomes more important on the background of weaker neutrino 
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Figure 18: Isotherms of the internal temperature Ti (or the surface temperature T^ , Table B) for a NS of the Geminga's 
age, as in Fig. M, but for the standard cooling (1.30 M©). 

emission produced by other neutrino reactions; the Cooper-pairing emission is asymmetric itself, being 
more efficient for neutrons than for protons. 

If the superfluidity is absent, a NS (1.30 M©) enters the photon coohng stage at t^ ~ 1.6 x 10^ yr. 
Thus the PSR0656 + 14 pulsar appears at the transition stage and Geminga is at the photon stage. 
The neutrino luminosity is already weak and the superfluidity affects the cooling mainly either through 
the heat capacity or through the neutrino emission due to Cooper pairing of neutrons and, to a less 
extent, of protons. 



Consider, for instance. Fig. 18. Two dashed auxiliary lines enclose the domain of joint superfluidity 
of nucleons. To the left of the upper line protons are superfluid and neutrons not, while below the lower 
line only neutrons are superfluid. The lines intersect at the isotherm of the temperature IgTj ~ 7.81, 
which a non-superfluid NS would have by the age t. The superfluidity with Ig Ten ^ 7.9 and Ig Tcp < 7.9 
does not appear by the moment t and does not affect the cooling. 

Horizontal segments of isotherms to the left of the upper dashed line show that the cooling is 
regulated by the proton superfluidity alone. For IgTcp > 7.9 K, the superfluidity appears just before 
the given moment t and is weak. In the range 7.9 < IgTcp < 7.9 + 0.5 it initiates the latent heat 
release and a weak increase of the total heat capacity (determined mainly by normal neutrons). In 
the range 7.9 + 0.1 < Tcp < 7.9 + 0.7 the weak raise of the heat capacity is compensated by the (also 
weak) increase of the neutrino emission due to pairing of protons. As a result, Tj initially increases, 
and then decreases. For IgTcp > 8.4, the proton heat capacity is reduced, and the total NS heat 
capacity decreases by ~ 25%. On the other hand, for high Tcp, the proton superfluidity onset is 
shifted to the neutrino cooling stage. The cooling delay produced by the suppression of the neutrino 
luminosity in the neutrino era is somewhat stronger than the cooling acceleration produced by the 
effect of superfluidity on the heat capacity. Thus Ti continues its growth with increasing Tcp. 
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Figure 19: Same as in Fig. |i|, but for a NS of the PSR0656 + 14 age. 

Horizontal segments of the isotherms end at temperatures T^n ~ 10^ K, at which the neutron 
superfluidity is switched on. The neutron pairing induces a splash of neutrino emission and cooling 
acceleration in a strip of the width ^^0.7 (in IgTcn) to the right of the auxiliary line. The minimum of 
Ti takes place in the interval 8.5 < lgTc„ < 8.7, because at 8.0 < IgT^n ^ 8.0+0.5 the neutrino energy 
losses are partly compensated by the latent heat release. The lowest temperatures Ti are realized in 
the case T^p S> 10^ K, in which the nucleon heat capacity suffers the strongest suppression. 

On vertical segments of the isotherms below the lower dashed line the cooling is regulated by the 
neutron superfluidity alone. With growing T^n in this domain, the temperature Tj varies in the same 
manner as with growing T^p in the domain of the purely proton superfluidity. This happens because 
at IgTcn > 7.9 the neutrino emission due to the neutron pairing is important. It speeds up the NS 
cooling in the range 7.9 < lgTc„ < 7.9 + 0.7 (see above). If 7.9 < lgTc„ < 7.9 + 0.5, the neutrino 
cooling is partly compensated by the latent heat release. Therefore the minimum of Tj takes place in 
the range 7.9 + 0.5 < lgTc„ < 7.9 + 0.7. It is not so deep as in the upper part of Fig. |l^ since the 
nucleon heat capacity is now suppressed only partly. 

A strong neutron superfluidity reduces the heat capacity stronger and the neutrino luminosity 
weaker than a strong proton superfluidity. Owing to the weakness of the neutrino energy losses this 
difference is sufficient for a NS with high T^n and normal p to cool in a different way than at equally 
high Tcp and normal n. The strong neutron superfluidity delays the cooling of those NSs which would 
be at the neutrino cooling stage or at the neutrino-photon transition stage if they were non-super fluid. 
This is demonstrated in Fig. |l^ (for PSR0656 + 14); in the absence of superfluidity at i = 10^ yr the 
pulsar would be at the transition stage. The strong neutron superfluidity accelerates the cooling of 
older NSs, e.g., of the Geminga's age (Fig. 18). 

Vertical segments of isotherms in Fig. 18 intersect the auxiliary line at temperatures T^p, at which 
the proton superfluidity appears. This superfluidity leads to the latent heat release and to the growth 



59 



of Ti in a strip of the width ~0.5 in IgTcp above the lower dashed hne. With further growth of Tcp the 
heat capacity is strongly reduced and the cooling is accelerated. For a very strong joint superfluidity 
of n and p {Ten ^ 10^'^ K, Tcp 3> 10^'^ K, the very right upper corner of Fig. 18), the nucleon heat 



capacity and the neutrino luminosity of the NS core are fully suppressed, and the cooling is governed 
by the electron heat capacity. 



Figure [19| is analogous to Fig. |18|, but corresponds to a younger NS. Its neutrino luminosity is 
somewhat higher and the relative contribution of nucleon-pairing neutrinos is smaller. The neutrino 
emission produced by pairing affects the cooling weaker. 

Comparison of the present calculations with observations will be made in the next section. 
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8 Thermal emission from neutron stars and superfluidity in their 
cores 

8.1 Thermal emission from neutron stars 

8.1.1 Observations of X-ray emission from neutron stars 

In the very first article devoted to NSs Baade and Zwicky (1934) predicted that NSs should be born hot 
and cool gradually emitting thermal radiation. Since NSs are compact, their radiation is weak. Modern 
detectors are able to register it only from the closest {D < 1-2 kpc) and sufficiently hot (Tg > 10^- 
10^ K) isolated NSs. The main radiation flux for the surface temperatures T<j mentioned above is 
emitted in the soft X-ray and hard ultraviolet spectral ranges (0.01-1 keV), which are thus most 
favourable for observations. This radiation, however, is strongly absorbed by the Earth atmosphere 
and can be detected only from balloons, rockets and space observatories. 

The attempts to discover NSs from their thermal surface emission undertaken in the middle of 1960s 
in the first balloon experiments with X-ray detectors were not successful. The NSs were discovered 
later, in 1967, as radio pulsars (Hewish et al., 1968). Soon afterwards they were also discovered 
in X-rays as X-ray pulsars. X-ray bursters and transients; however, radiation from these objects is 
determined by accretion of matter on NSs in binary systems, but not by outflow of thermal energy from 
the interiors of cooling NSs. A search for intrinsic thermal radiation from NSs was being continued. 

Wolff et al. (1975) as well as Toor and Seward (1977) attempted to detect thermal radiation of the 
Crab pulsar in a balloon experiment during lunar occultation of the pulsar. The thermal radiation 
was not detected but the upper limit of the surface temperature of the pulsar was established which 
does not differ considerably from the present upper limit. 

Subsequent observations were mainly carried out in the soft X-rays with the X-ray telescopes on 
board of the orbital observatories Einstein (1978-1981), EXOSAT (The European X-Ray Observatory 
Satellite, 1983-1986), ROSAT (The Rontgen Satellite, 1990-1998), ASCA (The Advanced Satellite 
for Cosmology and Astrophysics, operates since 1993) and RXTE (The Rossi X-ray Timing Explorer, 
launched in December, 1995). In these studies, various types of detectors have been used: gaseous 
(scintillation) photon proportional counters which possess moderate (from ~ 8 to 40%, depending on 
detector scheme and photon energy) spectral resolution and rather large angular resolution (~ 12"-3'); 
solid-state microchannel photodetectors with high angular resolution (down to 3") but giving almost 
no spectral information; solid-state CCD matrices with higher spectral resolution (to 3%), than the 
gaseous detectors, but worse angular resolution (1'), than the microchannel photodetectors. Temporal 
resolution (0.001-0.5 ms), in most cases, enables one to study variations of radiation fluxes concerned, 
for instance, with NS rotation. By the moment of completing this review the best angular, spectral and 
temporal resolutions possessed, respectively: the microchannel X-ray detector HRI (High Resolution 
Imager) on board of ROSAT (~ 5" in the energy range 0.1-2.5 keV), the detector SIS (Solid-state 
Imaging Spectrometer) on board of ASCA (0.4-10 keV, /S.E/E ~ 2% for photon energy £^ ~ 6 keV) 
and PC A (Proportional Counter Array) on RXTE (~ 1 ^s in 2-60 keV range). X-ray detectors have 
high sensitivity. For instance, the ROSAT observatory was able to detect the sources with a flux 
~ 5 X 10^^^ ergcm"^ s~^. More detailed information on active and future detectors can be found, for 
instance, through the Internet site http://heasarc.gsfc.nasa.gov/docs/heasarc/missions.htmL 

The orbital observatories have led to a great progress in search for the thermal NS radiation. Soft 
X-ray radiation has been detected from about thirty radio pulsars and from several isolated radio 
silent NSs. Thermal component has been identified sufficiently reliably in ~ 20% of cases. 

A search for the thermal radiation from the NS surface was also carried out with the ultraviolet 
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observatory EUVE (The Extreme Ultra- Violet Explorer, operates since 1992) using the DS (Deep 
Survey) and ST (Scanning Telescope) telescopes, sensitive in the photon energy range 70-760 eV, 
with the angular resolution ~ 30-45" (using wideband filters 40-190 A and 160-238 A). About 20 
known NSs were observed (Korpela and Bowyer, 1996). Five sources were detected in the 40-190 A 
band. Two of them, PSR 0656+14 and Geminga, are the cooling pulsars, and one, RX J1856-37, is 
an isolated (probably cooling) radio silent NS. For other objects, the upper limits to the ultraviolet 
fluxes were obtained. 

Since 1978 wideband photometric observations of isolated NSs have been carried out in optics 
using the sensitive ground-based telescopes such as ESO-NTT (The New Technology Telescope of 
the European South Observatory), the Keck Telescope, the 6- meter Telescope, and the orbital HST 
(The Hubble Space Telescope, launched in 1990). About ten radio pulsars have been identified as 
optical sources (Mignani, 1998). As a rule, the luminosity of isolated NSs in optics is several orders of 
magnitude weaker than in X-rays. This serves as an additional argument in identifying point-like X- 
ray sources as isolated NSs, especially if they are not pulsars. The point-like X-ray sources associated 
with accreting NSs in close binaries in pairs with normal stars (X-ray pulsars, bursters, transients) 
show usually higher optical fluxes due to contribution of companions or accretion disks. Multicolour 
photometric observations have been carried out for several isolated NSs (PSR 0540-69, PSR 1509-58, 
PSR 0656+14, Crab, Vela, Geminga), while for others (PSR 1929+10, PSR 0950+08, PSR 1055-52, 
RX J185635-3854, RX J0720. 4-3125) the fluxes have been measured in one or two filters only. Recently 
Martin et al. (1998), using the 10-meter Keck Telescope, obtained the first optical spectra (with ~ 2 A 
resolution), of Geminga, the isolated middle-aged pulsar. Earlier the spectroscopic observations were 
carried out only of two young pulsars, Crab and Vela, whose optical emission is non-thermal. Recently 
Shearer et al. (1998), using the MAMA (Multiple Anode Microchannel-Plate Array) detector with high 
temporal resolution and the 6-meter Telescope, discovered temporal variations of the optical emission 
(in the B-filter) from Geminga and PSR 0656+14, cooling middle-aged NSs; radiation pulsates with 
the NS spin periods. Earlier, pulsating optical emission has been observed only from several young 
pulsars (including Crab, Vela, PSR 0540-69) with the spectra of clearly non-thermal origin. The 
energy of the non-thermal emission is taken mainly from the NS spin energy. 

The optical and ultraviolet observations complement the X-ray spectra with the data in Rayleigh- 
Jeans part of the spectrum. This enables one to get more definite conclusions on the spectrum type 
(thermal or non-thermal) and on parameters of the NS atmospheres. Interpretation of the thermal 
radiation has to be consistent with observations of the non-thermal emission from pulsars and nearby 
space in radio, X-rays and gamma-rays. Finally, one should use modern atmospheric models of NS 
radiation for a correct interpretation of observations. 

8.1.2 Interpretation of observations 

Discovery and interpretation of the thermal radiation from cooling NSs is a complicated problem. As 
mentioned above, the radiation is weak since the NSs are compact. The radiation has to be separated 
from the surrounding background which is especially difficult for young NSs in supernova remnants. 
Strong background is created by non-thermal emission of synchrotron nebulae (plerions) which are 
formed around young and active pulsars. The background is produced also in non-thermal processes 
of emission of X-ray and optical quanta in the magnetospheres of radio pulsars. The distances to NSs 
are known poorly (uncertainty is often about some ten percent and higher). The next difficulty is 
provided by interstellar gas which absorbs soft part of the spectrum (from a fraction of keV down to 
optics), especially for the objects which are more distant than some hundred pc or are obscured by local 
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clumps of interstellar gas. Finally, the thermal radiation from the main part of the pulsar's surface 
is superimposed by the thermal radiation from the hotter (~ (2-3) x 10® K) polar caps heated by 
fluxes of energetic particles from the magnetosphere. The spectra of the both thermal components and 
the power-law spectrum of the possible non-thermal component strongly overlap which complicates 
interpretation even more. 

It is clear that the middle-aged radio pulsars (~ 10^-10® yr) are most favourable for studies of the 
thermal radiation. Their surface is still rather hot, Tg ~ 10^-10® K, to be observable but the objects 
are already insufficiently active to support synchrotron nebulae around them. Supernova remnants 
are dissolved by given age, and the non-thermal magnetospheric processes weaken due to pulsar spin 
down. All these factors reduce the background and increase the chances to discover radiation from 
the pulsar surface. The isolated NSs which, due to some reason, do not show pulsar activity from the 
early stage, for instance, due to the weak magnetic field or rapid spin down are of special interest for 
the cooling theory. These stars should possess lower background of the non-thermal emission which 
simplifies detection of the thermal component (see below). 

An observer which studies X-ray radiation from cooling NSs has usually a poor set of observational 
data in the form of photon count rates in various spectral channels of a detector; the count rates 
do not exceed, for instance, 1-10 counts per second for the brightest objects detected with PSPC 
(Position Sensitive Proportional Counter) on board of ROSAT. For the majority of isolated NSs 
observed with ROSAT, the count rate is ~ 10~^-10~^ counts per second; therefore, the exposure 
time should exceed ten hours to accumulate statistically significant number of counts (e.g., 100-1000). 
Data processing and reconstructing the real spectrum of the NS radiation is an example of a classical 
incorrect inverse problem which can be solved only under additional assumptions on the spectrum 
of the source. The spectral analysis is usually carried out by calculating the count rates in different 
channels by "transmitting" a model spectral flux of the object through a response matrix of a detector; 
calculated values are compared then with observations. A model flux depends on NS parameters such 
as distance, column density of interstellar hydrogen, etc. These parameters are varied (for instance, by 
X^ method) to obtain the best agreement with observations. However the response matrix is usually 
known with some uncertainty. Its parameters vary in time and require permanent checking using 
calibration sources on board of a space observatory or cosmic standard sources. The model spectrum 
may be multicomponent; for instance, it may contain contributions of thermal emission from the main 
part of the pulsar surface and from the hot polar caps as well as non-thermal magnetospheric emission. 
While fitting the spectrum one can, in principle, determine all the varying quantities: parameters of the 
NS thermal radiation (effective temperature, magnetic field and chemical composition of the surface), 
parameters of the hard radiation component (temperature and sizes of polar caps), spectral index and 
intensity of the non-thermal component, spin period, NS radius and mass, distance, column density 
of interstellar gas, etc. Additional constraints on the variable quantities can be obtained from similar 
fitting of the X-ray light curves. 

A complete realization of the above scheme is not possible, to date, mainly because of poor statis- 
tics of photon counts and poor spectral resolution of the detectors. Another obstacle comes from 
incompleteness of the theory of magnetospheric pulsar emission and the theory of atmospheres of iso- 
lated NSs with strong magnetic fields. The non-thermal magnetospheric radiation is usually described 
by a power-law spectrum with an index determined from fitting. In the absence of the atmosphere 
models, the thermal NS radiation is often described by the black body spectrum. However, properties 
of radiation emerging from a NS are actually determined by a thin atmosphere with the temperature 
growing inside the star. The atmospheric spectrum may noticeably differ from the black body. This 
is proved by the well elaborated models of NS atmospheres with "weak" magnetic fields B < 10^-10"'^'' 
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G (Zavlin et al., 1996; Rajagopal and Romani, 1996; Pavlov and Zavlin, 1998) and hot atmospheres 
{Ts > W^ K) with strong fields B > 10^^ G (Pavlov et al., 1995; Pavlov and Zavlin, 1998). In par- 
ticular, the radiation spectra from the weakly magnetized NSs covered with the atmospheres of light 
elements (hydrogen, helium) appear to be harder for the same effective temperature Tg. For typical 
temperatures Tg ~ 10^-10^ K, the light elements are almost completely ionized; the main contribution 
into opacity comes from free-free transitions which is well described by the Kramers formula (the spec- 
tral opacity decreases as E~^ with the growth of photon energy E). This makes the atmosphere more 
transparent for higher-energy photons, which emerge thus from deeper and hotter layers, leading to 
hardening of the spectrum. In the atmospheres of heavier elements (iron), which are more difficult to 
ionize, the photoeffect and bound-bound transitions are most important. This leads to weaker energy 
dependence of the opacity on photon energy. The radiation spectrum becomes closer to the black 
body, with the temperature corresponding to some mean depth from which the quanta emerge. In 
this case almost Planckian continuum contains deep and wide photoionization jumps and absorption 
lines of atoms and ions of heavy elements. In contrast to the black body emission, the atmospheric 
radiation is anisotropic even for the NSs with weak magnetic fields. The anisotropy depends on photon 
energy and chemical composition; the well known "limb darkening" effect is quite pronounced. The 
anisotropy is greatly amplified by a strong magnetic field. In this case, the anisotropy axis is aligned 
along the magnetic field and radiation beaming can be strongly asymmetric. An allowance for the 
anisotropy and gravitational bending of photon trajectories near the NS surface (Zavlin et al., 1995) 
is especially important for interpretation of the light curves. 

To date, among the models of NS atmospheres with strong magnetic fields, B > 10^*^-10^^ G, only 
the models of hot atmospheres (Tg > 10^ K) composed of fully ionized plasma are reliable. According 
to these models, the radiation spectrum becomes softer, than for weaker magnetic fields, but remains 
noticeably harder than the Planck spectrum. For lower temperatures, the magnetized atmospheres 
become partly ionized. The strong magnetic fields drastically distort structure of atoms and ions 
and affect the spectral opacity. Detailed calculation of the opacity is complicated and has not been 
done yet even for the simplest case of hydrogen plasma. The magnetic field increases the binding 
energies of electrons in atoms. For instance, the field B ~ 10^^ G amplifies the binding energy of 
the hydrogen atom to ~ 150 eV. The number density of neutral atoms becomes larger increasing the 
role of photoeffect and bound-bound transitions; these processes must be calculated with allowance 
for thermal motion of atoms. Such calculations are complicated since the electric fields, induced 
in an atom-comoving reference frame due to atomic motion across the magnetic field, are so strong 
that they change the structure of atomic energy levels and probabilities of radiative transitions (the 
dynamic Stark-effect; see, e.g., Potekhin, 1994). While calculating the opacity of this, partly ionized 
plasma it is important to take into account the effects of plasma non-ideality and pressure ionization. 
These calculations are being performed at present. Only the models of hydrogen and iron atmospheres 
neglecting the dynamic Stark effect have been constructed so far (Shibanov et al, 1993; Rajagopal et 
al., 1997). The dependence of the radiation spectrum on chemical composition has the same tendencies 
as for the atmospheres with weak field. The spectra of the iron atmospheres are softer and closer to 
the black body, than the spectra of the hydrogen atmospheres. An account for the induced electric 
fields, which strongly broaden the photoionization jumps and other spectral features, will enable one 
to check the above statement. 

Owing to the effects mentioned above an "atmospheric" temperature Tg appears to be noticeably 
(typically, by a factor of 1.5 - 3) lower than a " black body" temperature in interpretation of the same 
observational data. If the validity conditions of the available atmosphere models are fulfilled, then the 
spectra of the thermal radiation and the NS light curves are better fitted by the atmosphere models 
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than by the black body. Parameters inferred from atmospheric interpretation (emitting surface area, 
distance to a NS, column density of the interstellar gas, etc.) are, as a rule, in better agreement with 
the data provided by independent observations in other spectral bands (see, e.g.. Page et al., 1996; 
Zavlin et al., 1998). The black body interpretation often gives less realistic parameter values and 
meets some difficulties in explaining all set of observational data. 

Since statistics of count rates is commonly poor and the number of adjustable parameters is large, 
the confidence ranges of these parameters are often too wide. To make them narrower some parameters 
are fixed (NS mass and radius, distance, etc.). The parameters become then more constrained. In 
a sense, this is an illusion if not confirmed by the results of independent observations (e.g., in other 
spectral bands). For instance, radio observations of non-thermal pulsar emission are very useful 
for interpretation of X-ray observations. They give additional constraints on distance from parallax 
measurements, on interstellar extinction from dispersion measure determination, on orientation of NS 
spin and magnetic axes from polarimetric measurements, on pulsar age from measurements of spin 
down rate, on stellar mass from measurements of orbital parameters of binary systems containing NSs, 
etc. Additional constraints can also be obtained from optical and gamma-ray observations. 

8.1.3 Observational results 

As a result of more than 30-year search for thermal radiation, X-ray radiation has been discovered 
or upper limits of the surface temperature have been obtained for some tens of isolated NSs (Becker 
and Triimper, 1997, 1999). However, so far the discovery of the thermal component can be regarded 
as definite only in several cases. These cases include four closest middle-aged radio pulsars (Vela, 
Geminga, PSR 0656-1-14 and PSR 1055-52) and three sufficiently young radio silent NSs in supernova 
remnants (IE 1207-52 in the remnant PKS 1209-51/52, RX J0002-F62 near CTB-1 and RX J0822- 
4300 in Puppis A). All seven sources are reliably identified; their nature as isolated NSs is confirmed 
by observations of periodic pulsations (excluding IE 1207-52), produced by stellar rotation, and/or 
large ratio of the X-ray to optic luminosities. These NSs are listed in Table y in accordance with their 
age. For each star, we give its characteristic age (determined either from NS spin down rate or from 
morphology of a supernova remnant), spin period, method of interpretation of thermal spectrum, the 



effective surface NS temperature T^ = TgJl — (TZg/TZ) (Sect. 7.2.3) determined by this method, and 
distance to the source. 

Among the "atmospheric" interpretations presented in Table H, those for IE 1207-52, RX J0002-I-62, 
RX J0822-43 and the Vela pulsar seem to be most reliable. The interpretation of their spectra with 
the hydrogen-helium atmosphere models has an additional advantage: the fits of observational data 
with iron atmosphere models or with the black body spectrum are of lower statistical significance. 
Geminga and PSR 0656-1-14 are older and cooler. An important contribution into opacity of their at- 
mospheres comes from the effects of motion of neutral and partly ionized atoms across strong magnetic 
field (see Potekhin and Pavlov, 1997; Bezchastnov et al., 1998; and also Sect. 8.1.2). The presented 
interpretations have been made using the simplified atmosphere models in which these effects have 
been neglected. The hydrogen composition of the Geminga's atmosphere is additionally confirmed by 
a possible discovery of the proton cyclotron line in its optical spectrum (Bignami et al., 1996; Martin 
et al., 1998). The "atmospheric" temperatures of Geminga and PSR 0656-1-14 are likely to be closer 
to reality than the "black body" temperatures but less reliable than for younger NSs. 

Comparison of the surface temperatures of the objects from Table |^ with the simulations of NS 
cooling will be discussed in Sect. 8.2. In this discussion we will assume that the "atmospheric" and 
"black body" interpretations of observations are equally acceptable and obtain constraints on critical 
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Table 3: Observational data 



Source 


Igt 

[yr] 


P 

[ms] 


Modef^ 


[K] 


b) 


Distance 

[kpc] 


Refs. ^^ 


RXJ0822-43 


3.57 


75 


H 


6.23+^^ 


95.5% 


1.9 - 2.5 


Z 1999b 


bb 


6 61+0-05 


95.5% 


2.5-3.5 


Z 1999b 


IE 1207-52 


3.85 


— 


H 


^■WTol 


90% 


1.6-3.3 


Z 1998 


bb 


6.49+n? 


90% 


11 - 13 


Z 1998 


RXJ0002+62 


3.95^) 


242^^) 


H 


6 n3+°°3 


95.5% 


2.7-3.5 


Z 1999a 


bb 


6.i8t^:}^ 


95.5% 


3.1 + 0.4 


Z 1999a 


PSR 0833-45 

(Vela) 


4.4e) 


89 


H 


5.90tnt 


90% 


0.4 + 0.12 


Pa 1996 


bb 


f, 94+0.03 
0-24_o.03 


— 


0.5 


1995 


PSR 0656+14 


5.00 


385 


H 


r 79+0.04 


— 


n 9O+0.06 


A 1993 


bb 


c- qfi+0.02 


90% 


0.76 


Po 1996 


PSR 0630+178 
(Geminga) 


5.53 


237 


H 


c 9C-+0.08 


90% 


0.008-0.022 


M 1994 


bb 


cr 7cr+0.05 
"3- '^-0.08 


90% 


0.16 


HW 1997 


PSR 1055-52 


5.73 


197 


bb 


c OO+0.03 
0-O»_o,o4 


— 


0.9 


1995 



'^' "H" observations are interpreted with a hydrogen atmosphere model, 

"bb" with black body spectrum. 
' Confidence level of temperature T^ (95.5% corresponds to 2a level, a 90% corresponds to 1.64a); 

dash means that the level is not indicated in cited references. 
^' The mean age taken according to Craig et al. (1997). 

' According to Hailey and Craig (1995). 
^) According to Lyne et al. (1996). 



/) 



Z - Zavlin et al.. Pa - Page et al., A - Anderson et al., Po - Possenti et al, 
O - Ogelman, M - Meyer et al., HW - Halpern & Wang. 



temperatures of nucleon superfluidity in the NS cores. Here we will outline some other observations. 
The upper limits of Tg for the majority of other objects are too high to be of interest for the 
NS cooling theory. Let us mention only the upper limit of the surface temperature of the young 
Crab pulsar. Numerous attempts to detect the thermal radiation from the pulsar have failed due 
to the powerful non-thermal radiation from the magnetosphere and supernova remnant in the wide 
spectral range. The most stringent upper limit T^ = 1.55 x 10^ K was obtained by Becker and 
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Aschenbach (1995) from the ROSAT observations. Taking into account the uncertainty in distance 
and interstellar absorption, this upper limit does not contradict to the standard cooling of the Crab 
pulsar (with the non-superfluid core) although for stiffer equations of state the agreement is rather 
questionable. Let us mention that the surface temperature of Crab can be estimated independently 
from observations of after-glitch relaxation. According to the theory, this relaxation depends of the 
internal NS temperature, and the internal temperature is related to the surface one. Such estimates 
are not very definite but give T^ ss 1.6 x 10^ K (Alpar et al., 1985), in excellent agreement with the 
above upper limit. 

There are about ten more potential candidates to the cooling isolated NSs among radio silent 
compact X-ray sources in addition to those listed in Table ^. New observations are required to 
confirm that these sources are the cooling NSs. Two most probable candidates are the point-like 
X-ray objects in the young supernova remnants: IE 1841-045 in the center of Kes 73 and IE 161348- 
5055 in the center of RCW 103. The both objects have close black body temperatures, Tg ~ 0.7 
and ~ 0.6 keV. Such temperatures are too high for the ordinary cooling NSs of supernova remnant 
ages, ~ 2000-4000 yr (Gotthelf and Vasisht 1997; Gotthelf et al., 1997). However, for IE 1841-045, 
the above age agrees with its dynamic age determined from spin down rate (from observations of 
coherent pulsations of X-ray emission discovered recently), P = 11.8 s and P = 4.73 x 10~^^ s s~^ 
(Vasisht and Gotthelf, 1997). The standard (for radio pulsars) estimate of the magnetic field from 
P and P reveals the enormous field strength B ~ 10^^ G. The optical and radio emission from the 
objects has not been discovered. If confirmed in future observations (in accord with preliminary 
estimates by Heyl and Hernquist, 1997a), these data can be explained by the thermal radiation from 
young cooling NSs with superstrong magnetic fields and with outer shells composed of light elements. 
The presence of the superstrong field and the shell increases the thermal conductivity of the outer 
layers and, as a consequence, rises the surface temperature at the early cooling stage. For the second 
source, IE 161348-5055, this hypothesis should be treated with great care due to recently discovered 
mysterious raise of the luminosity for several months accompanied by the increase of the X-ray flux by 
one order of magnitude with the unchanged spectrum (Gotthelf et al., 1999). A search for pulsations 
with NS spin period has not been successful so far for this object, which has been observed regularly 
in X-rays since 1979. 

The sources of repeating soft gamma-ray emission (SGRs - soft gamma repeaters) belong possibly 
to the same family of cooling NSs with the superstrong magnetic field (called "magnetars" , Thompson 
and Duncan, 1996). However, they probably undergo different evolutionary stage. All of them are 
situated near the supernova remnants which indicates they are young. For two of them, SGR 1900-1-14 
and SGR 1806-20, coherent pulsations of X-ray radiation with the periods 5.16 and 7.47 s, respectively, 
have been discovered recently (Kouveliotou et al., 1998a, b; Hurley et al., 1999). Both sources have 
about the same large braking index as the source in Kes 73. Estimated surface temperatures of the 
objects {Ts < 1 keV) exceed typical temperatures of ordinary cooling NSs of their ages (t < 1000- 
10000 yr), if the ages are the same as for parental supernova remnants. 

Let us mention one more type of the sources. They are the so called "anomalous X-ray pulsars", 
the X-ray sources with the soft thermal spectrum (T < 1 keV; e.g., Mereghetti et al., 1998) and 
almost sinusoidal light curve. Typical periods of their brightness variation are about ten seconds. The 
magnitude of pulsating component of their radiation is anomalously small and more typical for the 
thermal radiation from the surface of a cooling NS than for an X-ray pulsar. One cannot exclude that 
at least some of these objects also belong to the family of cooling magnetars. No binary companions 
were found for these objects, and no accretion disks which would provide accreted matter for X- 
ray pulsars. A powerful magnetodipole radiation of a rotating NS with superstrong magnetic field 
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would explain naturally rapid spin down; an increase of the thermal conductivity of the outer layers 
would produce sufficiently high surface temperature. Young ages of some anomalous pulsars are also 
confirmed by coincidence of their positions with the centers of supernova remnants. They are also 
radio silent. 

Finally, let us mention two more radio silent isolated NSs. They have been discovered quite 
recently; first, as bright point -like sources of the soft X-ray radiation in the i? 05^4 T surveys, and later 
as weak optical objects in special observations with the ground-based Keck Telescope and the Hubble 
Space Telescope. These stars are not associated with the supernova remnants. 

One of them, RX J1856. 5-3754 (Walter et al., 1996; Walter and Matthews, 1997), has the black 
body surface temperature (8.5ibl.4) x 10^ K (Pavlov et al., 1996), for the adopted values ofTZ= 10 km 
and M = 1.4M0, and is situated at a distance D < 130 pc. Radiation fiuxes in the optical filters 
F606W (6000 A) and F300W (3000 A), detected with the Hubble Space Telescope, are 25.9 and 23.48 
stellar magnitudes, respectively (Walter and Matthews, 1997). They agree with the hypothesis on 
the thermal origin of the spectrum in the spectral range from optics to X-rays. The upper limit 
on the pulsating component is rather high (about 15%). Further accumulation of photon statistics 
will enable one to continue search for periodic pulsations. On the other hand, the observed thermal 
radiation can be explained by accretion of interstellar gas on a cold, weakly magnetized NS which 
moves slowly through interstellar matter (Walter and Matthews, 1997; Zane et al., 1996). In order to 
choose between these two hypotheses one needs to measure parallax and proper motion of the object. 

The other NS, RX J0720. 4-3125, is somewhat hotter (the black body surface temperature is 79 it 
4 eV). Its X-ray fiux is pulsating with the period 8.39 s and modulation ~ 12% (Haberl et al., 1997). 
There are not very reliable indications of the spin down, P ~ 10~^^-10~^^ s s~^. Optical observations 
with the Keck Telescope indicate the presence of an object identified with this source. The detected 
fluxes are 26.6 it 0.2 and 26.9 it 0.3 stellar magnitudes in the B and R filters, respectively (Kulkarni 
and van Kerkwijk, 1998). For a distance to the source D < 100 pc, the weakness of the optical 
emission is in favour of the hypothesis that the object is an isolated NS. One can expect to deal either 
with a cooling NS-magnetar or with an old, weakly magnetized NS accreting matter from interstellar 
medium. A choice of correct hypothesis can be done only on the basis of future observations. 

Let us mention also two more radio silent NSs, RX J0806. 4-4123 (Haberl et al., 1998) and 
IRX J130848.6 -^212708 (Schwope et al., 1999), discovered most recently in the ROSAT surveys; 
they are very similar to the two latter NSs. They are possibly the cooling NSs (Neuhauser and 
Triimper, 1999) but the observations are still insufficient to drew final conclusions. 

8.2 Theory and observations of cooling neutron stars 

8.2.1 Necessity of superfluidity for interpretation of observations 

The aim of this section is to confront the cooling theory of superfluid NSs (Sect. 7.2) with observations 
(Sect. 8.1). Table ^ presents the data on the surface temperatures of seven isolated NSs. In Fig. 2C these 



data (circles) are compared with the cooling curves (Sect. 7.2). Diagonal and horizontal shadings show 
the regions of the surface temperature T^[t) filled, respectively, by different standard and enhanced 
cooling curves calculated by varying Ten and T^p from 10^ to 10^'^ K. The dotted curves show cooling 
of non-superfluid stars. 

The "non-superfluid" curves are seen to be in poor agreement with observations. On the other 
hand, the observations can be explained by assuming superfluidity in the NS cores. This is illustrated 
by the standard cooling curve (the solid line) . The values of T^n and T^p are chosen in such a way to 
hit the maximum number of observational points at once. 
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Figure 20: Observational data on surface temperatures of NSs obtained (Table H) in the models of black body radiation 
(filled circles) and hydrogen atmosphere (open circles). Shaded regions show the ranges of T^ filled by the standard 
(diagonal shading) and the enhanced (horizontal shading) cooling curves of NSs with different Ten and Tcp (from lO'' to 
10^" K). The solid line shows the standard cooling (1.30 M©) for specified values IgTcn and IgTcp (given near the curve). 
The dotted lines exhibit the standard (upper line) and enhanced (lower line) cooling of a non-superfluid NS. The dashed 
and dot-and-dashed lines show, respectively, the standard and enhanced cooling of a NS which possesses an envelope of 
mass 7 x 10~^^Mq composed of light elements (see Sect. 8.2.1). 



According to Fig. 2C, all the "atmospheric" temperatures as well as the "black body" temperatures 
of RX J0002+62, PSR 0656+14 and Geminga are located in the allowed regions of the standard 
and enhanced cooling of superfluid NSs. Thus our cooling calculations can be compared with the 
"atmospheric" and "black body" temperatures Tg of these sources. 

High black body surface temperatures of RX J0822-43, IE 1207-52 (not presented in Fig. |2^, but 
given in Table y), Vela, and PSR 1055-52 are not explained by our models but can be explained by other 
models of cooling NSs, particularly, the models with superfluid cores. For instance, the observations 
of the Vela pulsar agree with the standard cooling of a superfluid NS {Ten = 10^ K, Tcp = W^^ K) 
possessing an outer envelope of mass ~ 10^^ Mq composed of light elements (Potekhin et al., 1997). 
High black body temperatures of RX J0822-43, IE 1207-52 and PSR 1055-52 may indicate either the 
presence of some additional reheating mechanism inside these sources (Sect. 8.1) or the presence of 
superstrong [B > 10^"^ G) magnetic fields (Heyl and Hernquist, 1997a). Finally, one cannot exclude 
that the black body interpretation of their spectra is incorrect (Sect. 8.1.2). 

8.2.2 Confronting calculations and observations 

Let us analyse which critical temperatures of nucleon superfluidity, Ten and Tcp, in the NS cores used 
in the cooling simulations (Sect. 7.2) agree with the observations of the NS thermal radiation (Table 
^. As seen from Fig. ^, we can explain the majority of NS observations either by the standard or by 
the enhanced NS cooling assuming the presence of nucleon superfluidity. These observations include 
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Figure 21: Confidence regions of Ten and Tcp which correspond to the "black body" surface temperatures (Table H, 
Fig. pol) of RX J0002+62, PSR 0656+14 and Geminga in the models of enhanced (left) and standard (right) NS cooling. 



six "atmospheric" interpretations (RX J0822-43, IE 1207-52, RX J0002+62, Vela, PSR 0656+14, 
Geminga) and three "black body" ones (RX J0002+62, PSR 0656+14, Geminga). 

Let us assume that the NS atmospheres may contain hght elements whose surface density is 
insufficient (Sect. 7.2) to affect the cooling. In this case, we can use the both, "black body" and 
"atmospheric" , interpretations of the spectra of thermal radiation. Although the values of M and TZ 
and, accordingly, the gravitational redshift for our models (Sect. 7.2.2) are somewhat different from 
those which have been obtained (or adopted) in interpretation of the observations, the temperatures 
T^ from Table |^ and Fig. 20 can be compared with our cooling curves due to several reasons. First, 
our calculations are not very sensitive to variation of a NS mass or radius (see Sect. 8.3 below). Second, 
the confidence regions of M and 7^, inferred from interpretation of observations, include, as a rule, 
the values we use in modelling. 

For simplicity, we assume that the internal structure of the indicated NSs is the same, and, partic- 
ularly, the NSs have the same mass. Then the critical temperatures of nucleons in their cores should 
be the same. Let us analyse the confidence regions of T^n and Tcp constrained by the observations. 
Including the observations of several NSs enables us to reduce these regions. 



The regions in question are plotted in Figs. 21-23. Figure 21 corresponds to the standard and 
enhanced cooling of NSs (with masses 1.30 Mq and 1.48 Mq, respectively) with the black body spec- 
trum. Figure 22 corresponds to the enhanced cooling of the stars (1.48 Mq) possessing hydrogen 
atmospheres. Finally, the standard cooling of NSs (I.3OM0) with the hydrogen atmospheres is shown 



In each figure, the lines of different types enclose the confidence regions of Ten and T^ 



cp 



in Fig. 

associated with the error bars of the observed NS surface temperatures T^ (Table ^). Correspondence 

of the lines to the selected NSs is displayed in Fig. 21. For the PSR 0656+14 and Geminga pulsars. 



the isotherms are taken from Figs. 16-19. In each figure the actual overall confidence region of Tc, 
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Figure 22: Same as in Fig. EU but for the enhanced cooling and "atmospheric" interpretation (Table 0) of spectra of 
thermal radiation from RX J0822-43, IE 1207-52, RX J0002+62, Vela, PSR 0656+14 and Geminga. On the right panel, 
we show in more detail the region (shaded rectangle) in which the allowed values of Ten and Tcp of five latter objects are 
either close or intersect. 

and Tcp lies at the overlap of the confidence regions of all objects. 



Figures 21-^ show that the observations of several NSs at once can be explained either by the 
standard or by the enhanced cooling, adopting either black body or atmospheric interpretations of 
the spectra. In all the cases, there are the ranges of T^n and Tcp close or joint for all NSs; they do not 
contradict to microscopic theories of nucleon superfluidity in NSs (Sect. 3.1). 

by adopting the enhanced cooling and the black body 



According to the left panel of Fig. 21 



interpretation of the observations we obtain two confidence regions of Ten and Tcp; each explains 
cooling of three objects at once. The first region corresponds to a moderate neutron superfiuidity 
(IgTcn ss 8.1) and a strong proton superfluidity (9.6 < IgTcp ss 9.8); the second, wider region 
corresponds to a strong neutron superfluidity (9.75 < lgTc„ ^ 10.0) and a moderately weak proton 
superfluidity (IgTcp < 8.3). For the standard cooling of a NS with the black body spectrum (the 
right panel of Fig. |2l|), there are also two regions of Ten and Tcp for the three NSs, but they are 
somewhat different than for the enhanced cooling. The first region corresponds to a moderate neutron 
superfluidity (8.0 < IgTc^ < 8.2) and a moderately strong proton superfluidity (8.7 < IgTcp < 9.1), 
while the second one is associated with a strong neutron superfluidity (IgTc^ > 9.3) and a moderately 
weak proton one (8.1 < IgTcp < 8.3). 



In the case of the enhanced cooling and the "atmospheric" surface temperatures (Fig. 22) there is 



the only region of joint or very close values of Ten and Tcp for flve NSs (the shaded rectangle on the right 
panel). It refers to a moderate neutron (8.2 < lgTc„ < 8.35) and strong proton (9.54 < IgTcp ^ 9.65) 
superfluidity. Such superfluidity (IgTcp = 8.3, IgTcp = 9.6) enables us to describe also the sixth object 
from the "atmospheric" set (RX J0822-43), if we assume that it possesses a thermally insulating 
envelope of light elements with the mass of 7 x 10^^'^ Mq (see Fig. |20| ). However, even in the absence 
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Figure 23: Same as in Fig. hi but for the standard cooling. 

of the envelope, the confidence region of Ten and Tcp for this object is sufficiently close to the joint 
confidence region indicated above. 

Finally, for the standard cooling and the "atmospheric" spectrum (Fig. |2^ ), there is again the 
only confidence region of Ten and Tcp, where the critical temperatures are nearly the same for the 
five NSs; it corresponds to a moderately strong superfluidity of n and p (8.25 < IgTcn < 8.35 and 
8.0 < IgTcp ss 8.4, shaded rectangle). The observations of RX J0822-43 can also be explained by the 
presence of the same superfluidity (lgTc„ = 8.3, IgTcp = 8.2) assuming that the object possesses an 
envelope of light elements of the 7 x 10~^^Mq mass (see Fig. ^). In the absence of the envelope, the 
confidence region of Ten and Tcp is outside of the joint confidence region since the error bar of T^ for 
this object lies at the boundary of the shaded region (filled by the cooling curves of superfluid NSs). 

8.3 Discussion 

As follows from the results of Sect. 8.2, the majority of NS observations can be explained using the 
adopted NS models by the standard and enhanced cooling for the "black body" and "atmospheric" 
interpretations of the spectra observed. In all these cases there are the values of Ten and Tcp close or 
common for all NSs which do not contradict to the microscopic theories (Sect. 3.1), but depend on 
cooling type and interpretation of the thermal NS radiation. 

The existence of the same critical temperatures for several NSs at once is quite unexpected. Initially 
we have expected different joint confidence regions of T^n and Tcp for different pairs of NSs. The result 
is even more surprising taken into account simplicity of our cooling models. It would be interesting to 
confirm (or reject) these results using more advanced cooling models and a larger number of objects. 

We have checked that our results are rather insensitive to variations of the NS mass M as long 
as M does not pass through the threshold value Mcr = 1.442 M© (of switching on the direct Urea). 
This fact is a consequence of the main simplification of our models, constancy of critical temperatures 
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throughout the NS core. In particular, the standard cooHng curves for the NSs with masses lower 
than Mcr practically coincide (see, e.g., Fig. |2|). The same takes place for the enhanced cooling curves 
of NSs (with masses M > Mcr), if protons or neutrons in their cores possess moderate or strong 
superfluidity with Tc > 10^ K (just as for the joint confidence regions of Ten and Tcp obtained above). 
We have assumed in our analysis that the NS ages t are known (Sect. 8.1.3). Then the confidence 
regions of Ten and Tcp for each star have been determined from an error bar of the surface temperature 
Tg obtained from the observations. Introduction of the error bars of At instead of the fixed values of t 
would lead to additional broadening of the confidence ranges of Ten and Tcp for each NS and to slight 
broadening of the joint confidence regions of Ten and Tcp. Let us mention that the characteristic ages 
of PSR 0656+14, Geminga and PSR1055-52 are determined within a factor of ~ 3. The uncertainty 
in ages of young NSs is higher which, however, should not cause stronger broadening of the Ten and 
Tcp regions due to weak slope of the cooling curves at t = 10^-10^ yr (see Fig. 2C). 



For the black body interpretation of the thermal radiation from NSs, our simulations predict 
a strong neutron superfluidity. Ten ~ 10^'^ - lO^'^''^ K, and weak proton superfluidity in the NS 
cores. This conclusion is in a qualitative agreement with the result by Page (1994). A strong proton 
superfluidity is also possible in the presence of a moderate neutron superfluidity. In such cases one 
superfluidity (n or p) is noticeably stronger than the other. In the both cases of "black body" and 
"atmospheric" interpretations of the observations the largest contrast of Ten and Tep takes place for 
the NS cooling with the allowed direct Urea process. Notice that such a superfluidity is predicted by 
Takatsuka and Tamagaki (1997b). 

In neither case we need simultaneously strong superfluidities of n and p to explain the observations. 
This indicates that the equation of state of the NS cores cannot be too soft (the softness would mean 
weak nucleon-nucleon repulsion at small separations which would induce especially strong pairing. 
Sect. 3.1). 

Finally let us mention that we can satisfy the observations by varying the only parameter, Tcp, 
if we assume the presence of the neutron superfluidity with Ten ~ lO^'^-lO^'^ K. For the standard 
cooling, this parameter should lie in the moderate range lO^'^-lO^''^ K, while for the enhanced cooling 
it should be larger, io9-45-io9-65 k. 
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9 Conclusions 

We have described one of a few methods to study fundamental properties of matter of density around 
the nuclear matter density or several times higher. The method is based on theoretical simulations of 
cooling of neutron stars and comparison of the results with the observations of the thermal radiation 
from isolated cooling neutron stars. 

We have described all the basic elements of the theory of cooling of neutron stars with the standard 
nuclear composition: stellar heat capacity (Sect. 3), neutrino cooling (Sects. 4~6), relationship between 
the internal and surface stellar temperatures (Sect. 7.2.1). However we have not intended to discuss 
all the cooling problems in detail but focused on a relatively new direction of the theory which had 
been initiated by the remarkable paper by Page and Applegate (1992). These authors were the first 
who noticed that the neutron star cooling could be very sensitive to the superfluidity of neutrons and 
protons in their cores. The temperatures of the superfluidity onset, in their turn, are sensitive to the 
equation of state of superdense matter. This opens the principal possibility to explore the nucleon 
superfluidity and, accordingly, the equation of state in the neutron star cores. 

The effects of superfluidity had been included into cooling simulations long before Page and Ap- 
plegate (1992). However they had not been the objects of special attention; systematic analysis of 
the effects of superfluidity on the neutrino luminosity and the heat capacity had been absent, and 
the effects had been described by the oversimplified expressions (Sect. 7.1). We have tried to fill in 
this gap and presented the results of detailed calculations of the heat capacity and neutrino emissiv- 
ities in superfluid matter in the frame of the Bardeen-Cooper-Schrieffer theory (Sects. 3 -6). Our 
consideration is done in a unified manner. The results are mainly presented in the form of simple fit 
expressions convenient for practical implications. Thus we have produced a rather elaborated method 
to include the effects of superfluidity into simulations of neutron-star cooling. This method can be 
extended easily to the neutron star models with superfluid hyperons. 

For illustration, we have presented (Sect. 7.2) the results of cooling simulations of neutron stars 
whose cores contain superfluid neutrons and protons. For simplicity, the critical temperatures Ten and 
Tcp have been assumed constant over the neutron star core. We have analyzed the effect of superfluidity 
on the enhanced and standard cooling and showed that this effect is really crucial. If the standard 
cooling of middle-aged (10^-10^ yr) neutron stars with non-superfluid cores goes much slower than 
the cooling enhanced by the direct Urea process, this is not so in the presence of superfluidity. As a 
whole, the nucleon superfluidity strongly reduces the difference between the standard and enhanced 
cooling. This enables one, in principle, to explain the majority of observational data by the standard 
cooling of the stars with superfluid cores. On the other hand, this circumstance disproves a widely 
accepted point of view that the direct Urea process necessarily induces rapid cooling. 

We have compared (Sect. 8.2) the results of cooling simulations with the observational data on 
the thermal radiation of neutron stars. Almost in all the cases (for different methods of inferring the 
surface temperatures from observational data) the observations can be explained in the models of the 
standard and enhanced cooling of a neutron star with the superfluid core. It is important that the 
required values of the critical temperatures Ten and Tcp do not contradict to the values predicted by 
various microscopic theories of nucleon superfluidity in superdense matter (Sect. 3.1). However our 
results should be regarded as preliminary; they cannot allow us to identify reliably superfluidity and 
equation of state in the neutron star cores. To impose more stringent constraints one needs to carry 
out more elaborated simulations of neutron star cooling taking into account density dependence of 
Ten and Tcp in the stellar core. It would be very desirable to make an equation of state of matter, 
used in cooling simulations, consistent with the critical temperatures. For this purpose, one needs a 
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representative set of equations of state of superdense matter (from soft to stiff ones), supplemented by 
the values of Tcn{p) and Tcp{p), calculated from the same microscopic models. This would enable one 
to construct a set of more realistic cooling curves and improve thus "calibration" of neutron stars as 
natural thermometers of superfluidity in their cores. Comparing theory with observations one would 
be able to constrain strongly the set of allowable equations of state. Unfortunately, simultaneous 
microscopic calculations of the equations of state and superfluid gaps are almost absent. On the other 
hand, the surface temperatures of the neutron stars cannot be determined uniquely from observations 
so far. 

We hope that after appearance of self-consistent calculations of the equations of state and superfluid 
gaps and after launch of orbital observatories of new generation {AXAF, XMM, ASTRO-E, SXG), 
which would produce more reliable data on the surface temperatures of neutron stars, this method 
will enable one to reach deeper understanding of the nature of superdense matter. 
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